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A BSTIiAG T 


The problem of designing an automobile transmission 
* 

an epicyclic gear box - is formulated and solved as an analy- 
tical problem and also as a constrained minimization 
mathematical problem* The design constraints that are 
imposed on the problem ensure that the number of teeth on 
any pinion are vithout undercut (or interference), the 
induced stress due to dynamic load on the teeth is less than 
the allov^able stress, when transmitting the maximum torque, 
the ratio of face width and circular pitch is within limits, 
the number of planets used ensuring assembly of gears in 
each train I all the gears in each train remaining coaxial# 
For the allowable maximum, torque, the overall size i.e. 'the 
volume is minimized, allowing a little variation in the 
desired speed reduction. , 

To show the usefulness of the proposed procedure, a 
comparison between the existing design of a transmission and 
that of the optimised one is made. 
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LIST OF. SYMBOT,;^ 


Suffixes 1, 2, 3 and 4 refer to the First, Second, 

Third and Reverse Trains respectively, 

Sj_,SNj_,NSj_ = No. of-teeth of sun gears', i = 1,2, 3, 4 

Aj_,Mj_,NAj_ = No. of teeth of Annuli, i = 1,2, 3, 4 

^i»^i?NPi = No, of teeth of planet gears, i = 1,2, 3, 4 

SKj_ = Ratio of facfe width to the circular pitch or 
Face width factor 

T = Maximum torque in kg-m 

FEND = Endurance stress of the material of the gears 

V,'\/ELSj^ — Pitch line velocity of sun gears in in/sec 
(or ft/min) 

hj_ = Face width of the gears, i = 1,2, 3, 4 
(XR)^ = Speed reductions, i = 1,2, 3, 4 

f = Static stress in kg/mm^ of gear materials 
fgt = Static stress in Ib/sq.in of gear materials 

= Total applied load per tooth in Ihs, i = 1,2,3, 
Wd^ = Total dynamic load in lbs, i = 1,2, 3, 4 
Wb^ = Bending load in lbs, due to endurance stress 
Cj_ = Deformation factor, i = 1,2, 3, 4 
ej_ = Error in action in Inches, i = 1,2, 3, 4 
Pj_ = Circular pitch of gear teeth, i = 1,2, 3, 4 
= Tooth form factor, i = 1,2, 3, 4 
W(z) = Objective function 
P(x) r: Penalty function 

[Hj Hessian Matrix or any positive definite matrix 



gi = Constraints, -1 =. 1 , 2 ,...., 25 ■ 

= G(I) 

* 

G¥(I) = Gradient of objective ' fmct ion 
GP(r) = Gradient of Penalty function 

R = Parameter weighing the penalty term of P- function 
s = Direction vector at the i'*^^ iteration 
- Step length ^ 

/j_ = Langrange 

(DIAS)^ = Diameter of sm gears, 1 = 1 , 2 , 3 , 4 
(DIAA)^ = Diameter of annulus gears, i * 1 , 2 , 3, 4 
= Module gears in the train, i = 152,3,4 
(PT)j_ *= No, of planets in tho gear train, 1 - 1,2, 3, 4 



CHAPTER I 


HTTROPnCTION 


The characteristics of a piston engine which is 
universally used in the passenger cars are such that it is 
necessary to make available a choice of reduction ratios, so 
that sufficient engine power may be obtained over a fairly 
wide range of road speed. The lowest ratio must be such that, 
starting from rest may be accomplished on the steepest 
gradients which are likely to be encountered; the higher ratio 
must provide for suitable acceleration to be maintained while 
speed increases and an appropriate top speed. These functions 
can be provided for by a manually operated gear box utilizing 
a conventional type layshaft gears or epicyclic gear trains. 

1.1 EPICYCLIC TRAIN 

An epicyclic train Fig. 1.1 in its simplest form . 
consists of three elements; a sun gear, an annulus (or ring 
gear) and a carrier, the later supporting two or more planet 
wheels which mesh with both the sun gear and the annulus. 

A difference in velocity ratio and/or direction of 
rotation may be achieved by holding stationary one of these 
elements and turning another, when the third member will 
respond relatively. The various speed ratios which can be 
obtained are described in Appendix B. 
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Epicyclic gears (planetaries) were widely used in 
the early days of motor cars and were introduced by F.W« 
Lanchester in about l899. These were superseded, with the 
notable exception of the Model T~Ford, by sliding layshaft 
gears (credited to Panhard in 1904). Again this type of gear 
was used in the well known Wilson pre-selective gear box. In 
1940, again a preference for the epicyclic trains is shown 
for the automatic transmissions. 

Epicyclic gears are particularly useful in automatic 
transmissions for the s3.mple reason that changes from one ratio 
to another can be made without closing the, throttle or losing 
traction, as there is no axial engagement or disengagement of 
gear teeth or dogs. 

Other advantages of epicyclic gears are as follows j 

1, An epicyclic gear train is extremely compact when compared 
with a layshaft train capable of transmitting the same torque- 
in fact, the length of an epicyclic gear train need be slightly 
more than the width of one gear face, the planet pinion bear- 
ings being contained within the pinion width. 

2, The tooth load in an epicyclic train is shared since there 
are at least two (usually three) equally spaced planets, 
whereas, in a layshaft train, the power is transmitted by a 
single pinion. 

3, A further advantage of having several planets is the 
balance, thus achieved, of the loadings on the bearings of the 



.3 

three elements, -whereas the shaft bearings of a layshaft 
train are heavily loaded. 

Thus an epicyclic train is inherently quieter since 
tooth loads are shared and, more Important, the pitch line 
velocity is reduced owing to the rotation of the carrier arm. 

4, The torque difference between input and output members of 
any gear train is ■ inversely proportional to the speed diffe- 
rence. In a layshaft gear train, if the output speed is, say, 
half the input speed, the output torque will be twice the 
input torque. The difference in torques is that reacted by 
the casing', which must be restrained from rotation. Now, to 
see the basic relationship between an epicyclic train and 
layshaft gearing, the reaction of the casing of the layshaft 
g^ar box will be the same as that on the carrier of the 
epicyclic train. 

5, A layshaft gear arrangement can, with a given set of 
parts, provide only two ratios,, but a given epicyclic train 
can provide six ratios * j two reducing ratios, two increasing 

ratios and two reversing ratios. 

' There are usually two limiting factors in the simple 
gear train design; one is the minimum diameter of the planets 
and the other is the minimum diameter of the sun wheel, which 
is usually imposed by the output or input shaft running through 
the centre of it. It can be easily shewn that it is impossible 


* Details given in Appendix B 
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to obtain a ratio of 2;i with a simple train. 

The compound trains are ideal for overcoming the 
simple planet ratio lijnitations and offers a great deal of 
freedom to the designer. By varying the ways nn which two 
gear sets are compounded and selecting input and output members, 
it is possible for the designer to obtain gear trains of equal 
ratio, but with high tooth load in one train and low in the 
other. 

For a given output speed, if the sun wheel is the 
output member, its speed will be low and so will transmit a 
high torque. Therefore the sun wheel should be avoided as the 
output member wherever possible as it results in higher tooth 
loads and a consequent increase in the size of the whole train. 

6. The gear ratio changes are made by applying hand brakes 
to cylindrical drums integral with the annulus. These occupy 
far less space as compared to sleeves, levers, etc. used in 

a lay shaft gear box. 

7. Another important advantage which epicyclic gears have 
over layshaft gears when considered for automatic transmi- 
ssions - the ability to change from one ratio to another 
without loss of torque transmission - is achieved by suitable 
clutches or brakes associated with the reaction members, one 
of which is released as another is engaged. 

8. Since a torque balance must exist when transmitting power, 
it will be apparent that the torque in the brake or clutch 
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holding a reaction member is the difference between input 
and output torque. Thus for large gear steps, the reaction 
torque is relatively high and it is a common practice to use 
a one-way clutch (or free wheel) or a self energizing brake 
band. For lower ratios, where the reaction torque is small, 
plate or cone friction clutches are more commonly employed. 

The potential speed of a reaction member is inversely propor- 
tional to the torque required to hold it stationary i.e, a 
reaction member which, by virtue of the gearing, would run 
at a high speed with the output member stalled, requires only 
a small torque to hold it stationary. The actual energy which 
is absorbed (and dissipated as heat) in a friction clutch or 
brake during a shift is related to the torque transmitted and 
the ratio step and is independent of the potential speed 
difference . 

1.2 LITBRATUBE SURVEY 

Quite a large number of papers on epicyclic gearing 
have been published. E.F. Obert^ describes the method of 
calculating input and output speed ratio from torque rela- 

Q 

tionships. L.A. C-raham gives the principles underlying 

3 

variable speed differential drives. A Verhoeef enlists 
advantages of some transmission mechanisms to be used in 
regulating transformers operated b-'’- remote control and for 
driving stirring devices in chemical experiments, when trans- 
mission ratio is high and pov/er to be transmitted is low. 
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C .Carmichael relates the -unlimited variety of possibilities 

in providing speed ratios bet-ween driving and driven members 

and simplification of calculation of speed ratios by routine 

procedure. V.Francis describes the possibility of using 

more than three planet pinions and formulas for spacing of 

pinions in the carrier. R»H .Macmillan^ gives the extra 

ordinary diversity of trains possible even -when only five or 

six -wheels are employed. A.Balogh^ gives the simplified 

formulas for calculation of velocity ratio of epicyclic trains* 
9 

M.A. Flint describes the importance of correct choice of 

numbers of teeth in sun wheel and annulus and equations for 

designing simple epicyclic trains with three pinions. 

10 

W.H, Maun gives a method for working out sirL5)le epicyclic 

11 

gear trains. S. Rapport explains the method of equal arcs 

for finding ratio of input to output in a planetary gear 

12 

train. D.B. Welboum gives a simplified method of calcula- 

13 

tion due to Beyer and Kutzback. R.N. Abild describes a 

simplified method for designing and selection of epicyclic 

gear systems, where choice of system is possible. E.I. 

14 

Radzimovsky expresses the method of finding efficiency, 

speed and power in following four advanced systems of planet 

gear, having complex kinematic relationships, two inputs 

15 

shafts or two output shafts. R.J. Willis explains a ■ 

method which permits rapid feasibility studies of epicyclic 

16 

gears. W.A. Tuplin enlists the practical possibility and 
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limitations of four-gear epicyclic trains in drive applica- 
tions and method of selecting the number of teeth of diffe- 
rent gears so that efficiency is m'aximum at high velocity 
ratios. R»M. Couhlin and James Ballmer-^ compute the forces 
and speeds for simple and compound gear trains, based on 
summation energy equations. H.G. Laughlin sfed A.R.Holowenho 
have developed the general design equations as functions of 
input power and expressions for functions are also establi- 
shed in terms of defined velocity ratio across entire system 

20 

and across control circuit. H.C. Toim relates the design 
of planetary gear and method of tabulating movements and 
analysis laws governing speed and direction of rotation. 

Most of these papers deal with the analysis of 
gearing of simple epicyclic trains. In no case any attempt 
is made to minimise the space occupied by coupled epicyclic 
trains which is characterised by two or more single trains, 
the rotating members of which are coupled together. Such a 
train is an indispensable sub-assembly in an automobile. It 
is certainly an advantage if the gear box designed occupies 
the least possible space, at the same time satisfying all 
the necessary and essential requirements. 



CHAPTER II 


DESIGN PROBLEM 

The system to be designed is a transmission gear box 
for an automobile . The function of this system is to deliver 
four forward speeds and one reverse speed to the output shaft 
(propeller shaft) from a single speed diesel engine for which 
the maximum torque is specified at a specific speed. The 
total volume of the gear box, excluding operational clearan- 
ces between adjacent gears is to be as small as possible. 

To show the effectiveness of the synthesis procedure 
a live problem is chosen. In one of the U.K. make buses, 'an 
epicyclic gear box is used. The gear box has the following 
specifications : 


Max. torque 62, i kg-m at 

1100 rpm. 

First speed reduction 

1 

: 4.28 

Second speed reduction 

1 

: 2.43 

Third speed reduction 

1 

: 1.59 

Top speed 

1 

: 1 

Reverse speed reduction 

1 

t 5.97 

Onr problem is to synthesise 

an epicyclic gear box 


of the above specifications in such a manner that of all the 
existing possible designs, the chosen one will occupy the 
least space. 

2.1 DESIGN CONCEPT 

There can be many schemes to satisfy the above 
specifications, but the scheme selected is that which is 
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actually TDeiug used in the exampla problem so that the synthe- 
sis made finds a good comparison. This scheme is shown in 
Fig, 2,1, The input is the engine side and the output is the 
propeller shaft. There are four epicyclic trains of gears so 
that three velocity reduction ratios, cno direct speed (these 
four being forward speeds) and one reverse reduction ratio, 
are obtained. The top gear drive is provided by directly 
connecting the input shaft to the output shaft by using a 
multiplate clutch. The first speed planet gear carrier acts 
as the driving member to the output shaft for all the forward 
speeds and the planet gear carrier of the reverse gear train 
acts as the driving member for the output shaft for the 
reverse gear. 

2.2 OPERATION OF GEAR BOX 

FIRST GEAR - The first gear speed is obtained by 
applying a brake to the first gear annulus so that it is held 
stationary. The engine is turning the main sungear so that 
the planet gears will be rolling round with it. This carrier 
is fixed to the output shaft through splines. Therefore, it 
imparts motion to the output shaft and so to the rear akle, 
SECOND GEAR - The second gear speed is obtained by 
holding the second train annulus stationary by applying the 
brake to it. The main sun gear being turned by the engine, 
causes the planet gears to revolve and turn their carrier. 

This carrier is connected to the first gear annulus which 
therefore turns, speeding up the rotation of its planet gears 
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and Carrier. In this case the output shaft runs faster than 
the first gear i.e. less reduction is achieved. 

THIBD GE5^ - The third gear speed is obtained by holding 
the third gear brake drum, which is integral with the sun gear 
running free on the engine shaft. The annulus of the third 
gear train is integral with the second gear planet carrier. 

The third gear train planet carrier is connected to the second 
gear annulus and drives it in the same direction as that of 
the engine shaft i.e._ increasing its speed. So the drive is 
taken back through the second gear planets and carrier and 
the first gear annulus, both of these are speeded up. This 
results in speeding up of the first gear planets and carrier 
which imparts faster motion to the propeller shaft. In other 
words by interconnecting the second and third planetary trains, 
an increase of speed is obtained at the first gear annulus, 
which increases the speed of planets and carrier* 

TOP GEAR - In top gear, all the trains are locked together 
so that they revolve as one solid cylinder driving the output 
shaft at the engine speed. This is brought about by releasing 
all the brakes and actuating the clutch which locks the sun 
gear of the third train to the driving shaft, through the 
drum. This causes all the sun gears to revolve at the same 
speed i.e. the engine shaft speed, 

RBYBRSB GEAk - The first gear annulus is connected to 
the sun gear for the reverse gear train and hence drives it 
in a direction opposite to that of the engine. When the 
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% 

brake is applied to the reverse gear annulusj the planets of 
this train carry with them the planet carrier in the opposite 
direction to the engine shaft rotation. As the planet gear 
carrier is directly connected to the output shaft, the direc- 
tion of rotation of the output shaft is reversed* 

2.3 DESIGN PARAMETERS 

It is convenient to consider three distinct types 
of parameters - independent variable, constant and dependent 
variable. The independent variables are those that the 
designer is free to select Independently in his design synthe- 
sis. Constants ai*e those quantities which are imposed by 
design requirements or have been decided upon in the first 
phase of the design process. Together, the independent 
variables and constants completely describe a design. 

Dependent variables are quantities introduced to 
. simplify the. mathematical relationships involved. Though the 
mathematical model can be constructed without the use of 
•these dependent variables, but they are needed to simplify 
the. mathematical relationships . 

.2*4 DESIGN SPACE 

If the number of independent design variables is n 
for a particular design problem, then every possible de.sign - 
can be thought of as being represented by a point in an 
n-dimensional space. This space is called design space... 
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2.5 CONSTRAINTS 

The constraints are limits imposed on individual 
parameters or on group of parameters, in order to insure 
their physical realizahility or their compatihility with the 
rest of the system and the environment. 

2.6 SYNTHESIS 

The synthesis of a design is the process which leads 
to the selection of a point In the design space. Any point 
in the space is called a ’design’ even if this gives impossible 
configurations. The analysis is carried out to know whether 
the design is worthwhile or not. 

2.7 OBJECTIVE FUNCTION 

If the designer desires to find the best design 
among all those in the feasible region, he must define a 
function f (x) , called the objective function whose value is 
a measure of the merit of the design. Conventionally f is 
designed such that the better the design, the lower its 
value. The best feasible design is the one for which fC x ) 
is minimum for 

g^ C X ) C for i ~ 1,2, m 
This formulation of the problem is called a mathematical 
programming problem. 
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COMPUTER AIDED DESIGTJ 
( DESIGN MALYSIS AND SYMHBSIS ) 

In this chapter, the important design constraints 
are discussed. For the salce of convenience in mathematical 
relations, all these constraints are expressed as inequality 
constraints. These are obviously functions of design 
parameters. 

The subscripts l, 2, 3 and 4 refer to first, second, 
third and reverse gear trains respectively. 

3.1 DESCRIPTION OF TEE PROBLEM 

For a given amount of energy to be transmitted it is 
quite obvious that the best system is that which is smallest 
in size. Decrease in dimensions leads to the decrease of 
weight, economy of material, reduction of cost, less vibra- 
tions, friction and wear. 

The volme is calculated by considering all the four 
gear trains placed side by side, as the space for operational 
clearances, carrier plates (arms), etc, is fixed and cannot 
be reduced. The volume therefore depends on the outside 
diameters and face width of the annuli of the gear trains# 

3.2 CONSTRAINTS 

The calculations are based on the following consi- 
derations or constraints : 




3.2.1 Condition of Coaxiality 


According to this condition, the axes of the central 
gears must coincide. This facilitates tapping energy from 
this mechanism or imparting energy to this. In Fig. 3.1, the 
arm A is the Inner and central gear 4 is the given element. 
Obviously in order that the gears 4 and 1 are coaxial 

^1 ^2 ^ ^4 “ ^3 C3*l) 

where r^, rgj r 3 and r^ are the radii of pitch circles of 
gears Ij 2, 3 and 4 respectively, 

3,2.2 Condition of Neighbourhood 

To minimize the load per planet and proper balancing 
of the arm, generally not one but a few planets are fixed to 
the same arm. The planets are generally placed at equal 
angular distances apart and they are placed on the same plane 
perpendicular to the axis of the gear 1 as shown in Fig. 3.2, 
The load per planet gear is minimum or the number of planet 
gears is maximum when their addendum circles just touch each - 
other as shown in gear 2 and 2*, From the triangle ABC, so 
that the addendum circles do not touch, the condition is 

2(r^+ rg) sin -^>2- rog 

where rOg is the radius of the addendum circle of the planet 
gears ‘and Q is the number of planet gears on the arm 

1T 
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rog = Tg + m, where m is the module 
T\ Z2 + S 


sin 


Zi + Zg 


(3.2) 


3,2.3 Condition of Assembly 

If Q is the number of equally spaced planets, the 
numbers of teeth in the concentric gears must satisfy the 
condition of assembly 
Z, + Zo 


Q 


= Integer 


(3.3) 


where z^ and Zo are tho numbers of teeth in the sun gear and 
1 

annulus gear respectively. 

The proof of this is as in Fig. 3,3, In this figure 

let P be the first planet to be assembled on the carrier R, 

the arrows indicate two teeth on P and A in engagement. If 

it is to be possible to assemble a second planet spaced at 
360 

Q degrees from the first, it follows that if R is rotated 
through ^ revolution with S fixed, the planet P must cause the 
gear A to rotate by an amount which brings another tooth into 
exactly the same position as that which originally engaged with 
P at its starting point. 

In ^ revolution of carrier R, the concentric gear 

Q 

A will move through 

Z3 


1 . 
Q 


. fs 

Z, 
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Z3 2l 

= QZ3 

For assembly of the next planet to be possible , this mst be 
a multiple of one pitch or ^ revolution of i.e* 


23 ^ 


^3 


= h = Integer, 


The induced stress in the material of the gears at 
the critical section of a tooth, should always be less than 
the allowable stress as explained below ; 

3.2'.4.1 Tangential Tooth .Load 

For a given horse power to be transmitted by a pair 
of gears, the tangential tooth load TiT is given by the follow- 


ing horse power equation : 

2 IT NT 


H.P. = 


4,500 


(3.4) 


where 


N - R.P.M. of the gear 
T = Torque in kg-m 


= WR 


W = Tangential tooth load in kg. 

R =r Radius of the pitch circle of the gear in m 
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3.2^4,2_. Beam Strengt h of Teeth 

The force ■which the tooth of one gear transmits to 
the other causes bending stresses in the teeth. The following 
Lewis equation gives the tangential load which the tooth will 
Carry in beam action : 

W = S^^bYP (3.5) 

where = Maximum bending stress in kg/cm^ 

b = Width of the tooth in the axial direction in cm 
Y = Lewis form factor which depends on the number of 
teeth in the gear and the system of gearing used 
= 0,154 - for 20® pressure angle. 

This form factor can also be found from tables on curves for 
a given number of teeth z. 

p = circular pitch in cm, 

3 ,2 ,4 .3 Dynamic Gear-Tooth Load 

The load applied to the gear tooth is greater than 

the transmitted load based on horse power. In general, the 

faster the gears are running, the more shock due to tooth 

errors and the more dynamic effects due to unbalance and 

torque variations in the driving or driven parts, 

k considerable amount of research has been carried out 

to determine the amount of dynamic gear tooth loads. A 

research committee of the American Society of Mechanical 

30 

Engineers - under the Chairmanship of B. Buckingham published 
the first authoritative work on dynamic loads. A simplified 
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formula giving a quick but approximate calculation of dynamic 
load, Wd is : 


Wd = W + 


0.05V(bc + W) 
0.05+ /be +W 


(3,7) 


where V = pitch line velocity in ft/min 
b =• face width in inches 
c = deformation factor - from Tables (30) 
or as given below : 

3, 2. 4,4 Determinatiop. of Deformation “Factor c 

It is necessary to have some measure of the error 
in action of the gears. The noise of operation is usually a 
very good test of the accuracy of gears. The following table^® 
gives some measure of the order of accuracy required at 
different pitch line velocities and should serve as a guide 
for the selection of the proper class of gear to meet speci- 
fied speed conditions. The values shown in the table should 
keep the noise of operation and the intensity of dynamic 
load within reasonable limits. 
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V ft/min 

T 

Error 

~ir 

iJL 

V ft/min : 

T" 

Jl. 

Error 

250 


p.0037 


2,250 


0.0013 

500 


0.0032' 


2,500 


0,0012 

750 


0.0028 


2,750 


o.oqi 

1,000 


0.0024 


3,000 


0,0007 

1,250 


0.0021 


3,500 


0.0007 

1,500 


0,0019 


4,000 


Q.0006 





4,500 


0,0006 

1,750 


0.0017 


5,000 


0.0005 





and 



2,000 


0,0015 


over 




Corresponding to these errors, the value of c is known from 
the table given below : 


VALUES OF c 


{Tooth § 'Error in a.ctlon~in inches ~ 

Gear Mate rial {Form 50.005 {O.OOl ro,002 50*003 5o.004 SO'.OOS 

1 - t, 5-- ... J I ....X-. : - L 


steel on 
Steel 


20® Fu 11 830 1,660 3,320 4,980 6,640 8,300 

depth 


3.2 ,4.5 Margin of Safety 

The static beam strength of the tooth should always 
be greater than the dynamic load as : 

For steady loads, = 1.25 Wd 
For pulsating loads ,W-j^=l,35 Wd 
For shock loads, W-j^^ 1*50 Wd 
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3l_ 

. Many gear designers have found that gears do not 
have to be designed to carry as heavy a dynamic load as given 
by Buckingham equation which gives the dynamic loada in the 
Order of l35 to 175' percent* 

Automative gears are usually fully hardened. Most 
of the time j they run at quite low torque loads. Frequently 
their operation at full torque amounts to lose than 100,000 
cycles. Tooth strength is usually more important considera- 
tion than wear. Dynamic overload due to tooth errors is 
usually not too important because tjie parts are light, the 
shafts are limber and the speeds are low. Misalignment 
effects are not too sepious because the face widths are 
narrow. Frequently the ends of the teeth ar® relieved on the 
tooth is crowned. This prevents serious trouble due to 
possible concentration of the lead at the ends of the teeth. 
So the above dynamic load type of calculations may be used 
as an alternative to the use of the velocity factor described 
below : 


3^2-,4»6 Dynamic Factor - Barth’s liquation 

The velocity factor is intended tq account for shock 
loading which results 'when accurate teeth are run together 
at high speed-. For spur gears the velocity factor or 
dynamic factor is given by t 


C 


V 


-C 


^ 01 -- 

3401+V 


(3.8) 
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where V - periferal velocity at pitch circle in in/sec« 

The allowable stress is then modified as given below in the 
Lewis equation : 

Sa = <=v (3-9) 

The Lewis equation is therefore written as 

W = Y b p (3.10) 

3,2.5 Ve.locitv Ratio Constraint 

The final velocity ratio can have a maxiimim varia- 
tion of 2 percent from the nominal values specified. This 
limit is based on experience. 

3jl 2...6 Interference Constraint 

To avoid interference or undercut in the teeth of 
the pinions, the minimum number of teeth should not be less 
than 14 but according to Buckingham^^ for epicyclic gears 
the minimum number of teeth should preferably be 16, i.e. 

16 (3,11) 

3.2.7 Planet Teeth 

The number of teeth, P, on the planets should be 

S 28 

equa.1 to and (A-^S) should he an even number • , 

3.2.8 Sun and Planet Centre Distance 

The sun and planets should be meshed at a centre 

QO ^ 

distance corresponding to (S t P + 2) teeth to avoid 
interference. 
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3,2,9 Annulll Inte rference 

The teeth of the annulii shoul'd be given a correc- 

' I 

tion coefficient of - 0.5 to avoid fillet interference when 
28 

assembling • 

3,2410 Tin Interference 

To avoid tip interference between the internal gear 

(annulus) and pinions (planets), the size of the pinions must 

oq 

be sufficiently smaller than the size of the internal gear 
3«2«il 'Face Width Restriction 

According to ISI ~ 3681-1966 the ratio of face 

width and circular pitch, k should be between 2 and 5. The 
best values considered from experience is 3 to 4. 

The module is defined as the ratio of the pitch 
diameter to the number of teeth. This must be always positive. 
Though the modules are standardised, these will be considered 
as continuous variables. It is usual, from experience to 
have the \alues of modules not greater than 10, The preferred 
values arc : 1, 1.5, 2, 2.5, 3, 3,5, 4, 4,5, 5, 5,5, 6,6.5,,,.. 


3,2,13 Size of Annulii 


(3,12) 


The diameter of the annulus of first, second, reverse 
and the drum for the third gear and top speed should not have 
much variation to get a compact gear box. So the variation on 
the diameter allowdd to the annulii of second and reverse 
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gear train is restricted to 5 percent with respect to the 
annulus of first gear train. As the third gear train annulus 
is enveloped by the second gear annulus, allowing free move- 
ment of each, so the variation allowed is such that the dia- 
meter of third gear annulus is 15 to 5 percent lower than the 
diameter of the first or second gear annulus i.e. 


1.05 DIAA1^,DIAA2 0,96 DIAAl 

(3.13) 

0.95 DIAAl ^. DIAAS ^^0,85 DIAAl 

(3.14) 

1.05 DIAAl ^DIAA4 ^ 0.95 DIAAl 

(3.15) 


3.2,14 Machining Cost Reduction 

To reduce the cost of manufacture, the sun gears of 
the first and second gear trains are made common in one 
integral unit which implies that the module will be the same 
for the first and second gear trains. 

3 .3 ANALYSIS FOR OF TEETH 


By considering the condition of coaxiality and the 
module of the gears in a train to be the same , for the given 
velocity ratio (XR)j_, the tabulation method is used for 
calculating the relations between the teeth of various geari^, 

^*3 ,ri ;r a. t ,, „'r j.al£t„(Fig.3.4) 

With annulus held stationary with braJce 


— — 

1 — ±1 

6 step 2 5 0 

1 « 

-4— - 

vResult a +1 

2 2 

2 L 



T 

t 

i 


4 

5 
2 
i 
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A (XR) -3: 


1 + ^ 


or A = , (XR) - 1 i s 
” 1 j 1 


. fiiii 


L 2 


3tg_Llf2 Second Train (Rlg.3^) 

Annulus is held stationary. 
For train 2 


(3.16) 


(3.17) 


c 

5 

T 

5 

5 Sp 

I 

{Step 1 

5 

0 +1 

X_ 

1 

1 

fi 

fi 

OStep 2 

5 

0 

0 0 

4 

2 

2 -1 
1 

5 ^ 

I " ^2 

“1 

5 

5 

fi 

5Re suit 

11 

T 

0 +1 
i- 

5 

0 0 

5 

8 sg 

5 

—A 


Now = 1 = Rg- 
\ = ^2 

i' * 

R^ revolves at x revolution 
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0 R 

4 - ^ 

0 A (=Ro5 

01 ^ 

0 s4=Sp) 

0 1 2 


4 

step 1 

11 

i X 

5 

0 X 

0 

8 

6 X 

^ 0 


5 

4 

Step 2 

{ 

0 0 

8 

0 l-x 

8 

0 -d-x) 

h 

,q 

1f 

5 


J 

_J 

5 


4 


5 Q 0 XTfi 

Result 5 2 0+1 0 x( g ■ 0 

0 5 d i-fil 


Equating = S 2 


Si ^ - Si 


Or 


= 1 ^b— 

! Ai+Si J Sg ^ Ai+^i 


A *^S 

■/. Gear Ratio (XR)„ = -—;— 

$2 •x 


(3,18) 


X = C22 




= So X CXR)o-S 


2 


= S„(C22) (XR)p-S. 
2 1 


3, 3, 1,3 Third Gear Train (Fig.3.6) 


S is held stationary 

Cj 


(3.19) 


0 

5 

|step 1 


T 


% 


+1 


S3 


+1 


5 

5 

T 

0 

0 




+1 


S 

0 


step 2 


0 

5 

JResiat 

5 


0 


+1 


-1 


0 


^3 


. ^3 0 

II ~^30 
^ 0 



For train 2 ■ Ag. = R 3 , Rg " 


jj r-- 

0 

il 

0 Step 
ii 
0 

0 Ptap 2 0 0 

6 .5 

0 0 

5 Result 5 G- 

5 -3 


T~ Rg (-A3) 


-4- 


i*g(=R3) 


G. 


I'-Gg 


+1 


0 

j 

If" 

] 

L. 


^2 


5 

0 

0 
5 

8 

Ao 8 
41 -G 3 ) ~ 8 
___r2 8 


G. 


'3 


GoC 


■‘'^2 "^2 


V^SgJ 


In train 1 speed be x-unlp.om 


8 0 

3 8 

i .. ... 0. 

R 

fi 

3 

A-(=Ro) 

X ^ 

8 

8 s (=Sp) 

“8 

8 

8 8 

3 step 1 0 

1. .,4 

X 

3 

8 

X 

8 

8 x 

8 

8 

->_8 

8 8 

8 step 2 8 

J— S 

0 

8 

8 

8 

G3-X 

8 in, 

8 •-(G3-X) ^ 

- 8 ^. 1 . 

8 

8 

_8 

1 i 

8 Re suit 0 

i J 


8"~ 

fi 

) 


0 A,+S, G„ 

5 

3 ®1 ‘ 

:l{ 


Sr 


‘■'3 


By equatt ; uie results oi’ ^nd S 


^3 ^ 


/ig + S 


2 


’2) > ^ „ 13^-1 


S. 


1 


Go A., 


^ . 0 ,zi_ ) I ti . is41?.] _ Mi — 
^ - ^3 •a^+s^4_ s; ^2 J Sgdi.+s^) 


C3.20) 


(3.21) 
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.. Gear Ratio (XR)g 







1. 


and 


t. 




•''•2+S2-, 
So ^ 


X 


S 3 

G 3-1 




(3.22) 

(3.23) 


... Top ..Qegj; (Fig -3. 7) 

'The top gear speed (l;l) j.c obtained by me.ans of 
the mnltiplate clutch. By engaging the multiplate clutch, 
the sun gear of the third train sun gets connected to the 
input shaft r. So the sun gears of first, seccnd. and third 
train run at the same speed simultaneously. So uhe whole 
system moves as a solid cylinder at the engine speed which 
is directly fed to the propeller shaft*, 

3 .3. 2 Rever se Sueed (Fig, 3, s'* 

E.’ake is applied to A_^ to hold it stationary 


J Step 2 


(Result 


R 4 

5 ft : 

0 fi 

0 1 5 

^ 5 

? C? 4- ^ O 


0 


5 

0 

.X„. 


1 


9 

— ft 

) 

_JL_- 

0 

0 

..J 


+1 


0 


( 

> » 

.L-_ 

0 

0 


-i-1 


A/ 


1-'- '^■4 
1 , . =w^ 

. 5 


5 

5 

0 

A 

S 

5 
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Taking = Sg 
C22 = 


+ 1 




_ S^[i-(XR) 2 j (XR) 
" TxRTg^'CXR)^ 


(iii) III Speed 


If ( 


ll. 


3 I^j+S^ 


'h , 

.^1 Sg . 




SgCv 


‘'■2*^2s h. 


CgCXRjg = Gg (-1^) . ^ 


which give 


or 


3 


(XR)3A^S2 + (XR)2S^(ii2+Sg) - ( A^+S^) (i^+Sg) 
So 


(3.29) 


and 

=3 

= 1 + 

or 

^*3 

_ ^3 

G 3 - 



r (XR) 

or = 




(3.30) 


ftl 
3 A2 


""2 


(XR).(I-'P-) - (XR) 


L 


It 


1+ (l+7r“) 




(XR)^-l) 


,.(3.31) 




(XR)^ , (I 


“^3 " 


3 1- 


.R) 3^-1] [(XR) 2- (XR) 3J + vXR; ^^[(XR) (XR) J|[(XR) 3_(XR)3-: 


CXR)3 "(XR)i |■l-(XR)2j-{(XR)^-lj{(XR)2-(XR)^’J~^SR)^ 




(3.32) 
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(iv) Reverse- Speed 
C 4 = I + 



or 





= 1 - (I-C4) ^ 

laE)^ -1 j S4 

(XE)j^ -1 


(3.33) 

(3.34) 


3.5 SUMMilRY OF CORSTRiin'ITS USED DI PRO GRiU^Ii^IB 
1 * Interference constraints 

2. Neighbourhood constraints 


• ^ Tt \ 


Sin 


3. Stress constraints or Horse Power (Torque) Constraint 
(1®)^ ^ 1.35 (WD)^ 


\ 


2T 


~T"F“ 

h m. Y_. 


1 1 


4, Size constraints 
DI/J.2^0.95 DliJll 
^1.05 DliJll 
DIMS ^ 0.85 DI/ull 
^0.95 Dl/Jil 
DIAA4^ 0,95 Dl/il.! 
^1.05 DIA/il 
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5. Assembly constraints 

+ Si . 

= Integer 

6. Velocity -Ratio Constraint 

1X2 (ZR). '> (XR)_. ^ 0X8 (ZR). 

7. Interference constraint for plrnet gears 

(FW). \ 16 

8. Module constraints 

^ 1.0 

3 . 6 COMPUTER PROGRAMMIN G 

The programming for the design has four different 
stages. The output of the first gear train together with the 
variables of the second gear train, is the input for the third 
gear train. Similarly the input of the reverse gear train 
consists of the output of the first gear train and the velocity 
ratio of tb ;■ reverse gear train.. The flow chart given below 
explains ti ■' process of l^'^ration in deciding the number of 
teeth of aT’ luIus and planet gears ^ the murber of planets, 
module and diameter of annulus. 



start 


INFDT :T ,FSm) ,RH4 ,XRl, Si 
For I Gear Train 


Round off At , Pn to Integer Values 


Commute Lewis Form r actor 



Compute Max .No .of Planet | 


_] Set No .of Planets j 


JT PN1< 2^ 
"''^^ 00 ^ 


If neighbourhood condition 
^ satisfied — 


f assembly condition satisfie 


Set Module fr^g pre^rred numbers 


Compute DLxSlj DIj,‘iA1 


If DIA/a 300.0 mm 
fYes"" ^ 


If bending load due to ^ 
— ^ -^1.35 Dynamicio 


PRINT SNl, INI, PKl, PTl, AMI, XRl, DIAAl 


t 


No 

soluti 



Continued 


















1 1 


Input : SHI ,xilTl,i:J‘dl ,DI1A1 ,XRl ,XR2 
FOR II GR.Y1 TRilF 

! 


Set XR2 = 2.41 + kl 


j Compu.te A 25 P 2 j 

^ 

Roimd of to into go rs ' 

. , -1 - L*fm -I , _ -- - ^ ^ 



Compute Louis Form Factor ~| 

- 

Compute maz.iio.of planets ■ 


[R2 Xfl2 
(ma3S 



Set Ho. of Planets = H'ILX2-BH 


:f BT2 2 


If neighbourhood condition 
s satisfied — 


If assembly condition 
_ satisfied 


Set module from preferred 

nos. iH2 = Lin | 

V 


Compute DIAS2, DI^J^ 

I 



L 


ho 

solution 



If 1.05 DI.Ja DIih.2 C.S5 DIiJ.l 

Yes 


If bonding load duo to FH''!) 1,35 

— D;''namic 


I Yes 

I PRINT SI'I2, A^I2, HI2, iii'i2, 2R2 , DIixA2 


Stop 


Similar flow diagram for III and reverse gear trains 
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3.7 RESULTS 

Starting from a particular number of teeth of sun 
gear of first gear train, we get three sets of possible com-* 
binations of teeth of second, third and reverse gear trains^ 
Minimum and mazimum number of teeth of sun gear of first 
gear train were restricted to 16 and 36 respectively. Actual 
number of possible solutions depends on the variation, of 
individual sizes. The actual number of solutions depending on 
the various tolerances on annulus diameter, is given in 
Table II. 

In view of the innumerous number of solutions, when 
tolerance zone is rather slack, the results forabou.t one percent 

variation on the annulus diameter of different gear trains 

■ 

ar^ presented in Table I. 

It is quite obvious that if one eleme]||; is picked 
from each set, a design will follow to obtain a 
design, which gives the minimum volume of the, gear box for 
any given number of teeth of the first sun gear. 

Tbe criteria of selection of the minimum size is 
the minimum of the function : 

4r 2 

(DIAA) (iM). 

1=1 i ^ 

(The face width being directly proportional to the module). 

The results are presented in Table I. 



TABIE 1 


si. 

flNo. of 

. PIRST 

mmn 


No. 

Jteeth/ 52311 

Oinodule of 5 
flsun of I 5 
kear train 5 

SNo.of ilHCllDIiU.lf 'bl .0XR2 

Heeth i_ . 5 5 fi 

Oof 5 5 0 5 

_5a15 si5 5 5 5 . . 

lT]~PTfMMi|l'i| l| II 1 

1 

21 

4.28 

69 24 

5 

225 

39.0 2.43 21 

'69 

24 3 

225 

2 

16 

3.5 

4.28 

62 18 

4 

182 

38,5 2,42 

16 

3.5 

52 

18 4 

182 

3- 

17 

No solution 








4 

..18 

- No solution 








•5 

• 19 

3*5 ■ 

.4.28 

62 22 

3 

217. 

38*6 2.42 

19 

3.5 

62 

22 3 

.217 

■6 

20 

-3*5 ' 

'.4.28 

66 23 

2 

231 

38.5 2.43 

20 

3,5 

66 

23 2 

231 .. 

7 

21 

3,0 

4,28. 

69 24 

5 

207 

33.0 2.43 

21 

3,0 

69 

24 3 

207 

8 

22 

3v5 

4.28 

72 25 

2 

252 

38.5 2.42 

22 

3.5 

72 

25 2 

252 

9 

23 

No solution 








10 

24 

No solution 








11 

25 

No solution 








IS 

26 

3.0 

4.28 

85 30 

3 

255 

33.0 2.42 

26 

3.0 

85 

30 3 

255 

13 

27 

No solution 








14 

28 

2.5 

4.28 • 

92 32 

5 

230 

27.5 2.42 

28 

2,5 

92 

32 5 

230. 

15 

16 

4,0 

4*28 

52 18 

4 

208 

44,0 2.42 

16 

4 

52 

18 4 

208 

16 

21 

3.5 

4.28 

69 24 

5 

241 

38.5 2.43 

21 

3.5 

69 

24 3 

241 
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1.59 24 

60 

18 

3 

200 

38,0 

5.97 

32 

68 

18 

5 

225 

39,0 

5960 

I 

1 

1.61 21 

55 

17 

4. 

165 

33.0 

5.97 

35 

75 

20 

5 

187.5 

27,5 

4850 

£ 


3.0 2.5 

No solution 


No solution 


1,57 

19 

3.5 

56 

19 

3 

196 

38.5 

5.97 33 
3.0 

71 

19 

4 

2l3 

33.0 

5550 

1.61 

23 

3,5 

61 

19 

3 

213,5 38,5 

5.97 36 
3,0 

76 

20 

4 

228 

33.0 

6370 

1.61 

24 

3.0 

64 

20 

4 

192 

33.0 

5.97 28 
3.5 

59 

16 

3 

206.5 38.5 

4630 

1,57 

22 

3.5 

65 

21 

3 

227 

38 #5 

6.01 39 

3 

84 

22 

3 

252 

33,0 

7400 


No solution 


No solution 


No solution 

1,57 20 58 19 3 232 44.0 5*96 30 64 17 4 256 44.0 7540 

4.0 4,0 


No solution 


1.58 

24 

3 

69 

23 

3 207 

33.0 

5.96 

36 

3.0 

76 

20 

4 

228 

33.0 

4940 J 

1.59 

23 

3.0 

64 

20 

3 192 

33.0 

5.97 

28 

3.5 

60 

16 

4 

2l0 

38,5 

5550 « 

1,59 

24 

3^5 

64 

20 

3.5 224 

38.5 

5,97 

32 

3.5 

68 

18 

5 

241 

38.5 
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TABLE II 


No. of teethfi 

r 

5 


POSSIBLE NBMNER OF SOLUTIONS 


in 1st Gear 


Percentage variation of Anniili- diameter of 2nd, 
Rpyejss. .I r vin s •w .frJTT.tt. I,g t Tya ift 


t 






16 

3328 

17 

- 

18 

- 

19 

252 

20 

208 

21 

1800 

22 

51 

23 

- 

24 

- 

25 

- 

26 

54 

27 

- 

28 

78 


22248 


7976 

2304 

8904 

2808 


1250 


78 


CHAPIER IV 


MATHEMATICAL PR0GRA^D4IHG APPROACH TO TOE p-FjQ-RT.T^ 

In this chapter} the sane problem Is solved as a 
mathematical programming problem. 

The general mathematical programming problem is to 
determine a vector x that minimizes the function f( x ) 
subject to the constraints (x)^0, i = 1,2 , ^m. For 
this design the vector x becomes a vector of design variables* 
The function f (x) , called the objective function, with 
respect to which the design is optimized as a computable 
function of the design variables, which expresses the quali- 
ties that make one design better than another. The constraints 
g^ (x).^0 are the design restrictions, the satisfaction of 
which distinguishes acceptable designs from unacceptable ones. 
The design ■which optimizes the objective function (i.e, minimize 
the volume) and satisfies all the design constraints is ’’the 
best design". Because of the non-linear characteristics of 
■the objective function anii the design constraints and because 
of the special nature of certain constraints, the method 
seleptpd is the unconstrained minimisj^tion method which uses 
penalty fuhhtions to handle the constraints. 

4,1 MATHEMATICAL PROGR^NTMTWn PBOBT.-RM 
4,1*^ The Mathematical Model 

Here the mathematical model of the gear box is 
formula^ted. The tejm mathematical model is considered here 
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to mean the complete set of equations and inequalities, , which 
are necessary for the statement of the design problem as a 
mathematical programming problem. This set of functions 
includes all of the equations and inequalities needed to des- 
cribe the kinematiCj, geometric and dynamic requirements, 
including coaxiality, neighbourhood, assembly, interference, 
d 3 niamic load, etc., the design function and design constraints. 
4 * 1 . 2 — The Objective Function 

The object of the synthesis problem being considered 
here is to design a gear box using epicyclic gear trains which 
transmits the given torque and gives the designed speed ratios 
within 2 percent of the specified values. The values pertain 
to a heavy truck/bus having four forward speeds and a reverse 
speed. 

These following twelve quantities are taken as design 
variables : 

Modules, m^^ j i = 1 j2,3,4 

Velocity Ratios (XR)^'=B^ , i = 5, 6, 7,8 

No, of teeth on Sungears S^^B^ , i = 9,l0,ll,l2 

Any set of values for these design variables is called a 
’design' even if it is absurd (i.e. gives gears which cannot 
be assembled) or inadequate in terms of transmission of spe- 
cified energy etc. 

The synthesis problem may now be stated more specifi- 


cally as follows : 
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An epicyclic gear box is to be designed so that 
1) it transmits the required torque 
ii) it gives four forward speeds and one reverse 
speed within a variation of 2 percent. 

The volume function can now be fbrmulated as the sum of 
volumes of the four gear trains, the top speed being obtained 
directly. The volume that a gear train occupies, depends on 
the outside diameter of the annulus and face width of the teeth. 
The space occupied by the pins, arms, etc. and multiplate clutch 
do not directly effect the volume function as these occupy 
a specific space and so are neglected in the optimization 
problem. ^ 

so f (55 = 5 f (4.11) 


The best design is the one which minimizes this total 
volume and satisfies all the design constraints. The design 
constraints are the restrictions placed on the minimum number 
of teeth, which can be cut without interference, the maximum 
number of clanets, the dynamic load tc remain less than the 
load calcilated from endurance stress (or the load calculated 
by using velocity factor), and the number of teeth on annulus 
and sun gears to affect the assembly. 


4,1,3 Objective Function Formulation 

The object of the optimization problem is to minimize 
the volume of the gear box, which is given by : 
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^ of cross section of JUinulus X face width 


'ff o 

“ ^ ^ ^i“^i » ^ ~ 1,2, 3, 4 


T]'"^ ^3 2* 

= ~ ^ nii X (4.12) 

The objective function can now be written as 

iT^ - 

^ = 4^ ^! { “l {“l (c22S2XR2-S1)^| 


where 


^(“s B s|] 


U) 

_ ,■‘'2*^2, 5a ■'h 
=22 = < 1^3 Sa 15^ 

S-1 r An 1 ApS-1 

foo - n^.( 4. .-A ■ .gf , — ,, V, 1 

A^+Si [ ^2 J 


(4.13) 


^3 =1^1; 

C44 = 1 + A4/S4 

4.1. 4 The Design C onstraints 

The design constraints or design restrictions, a 
satisfaction of which, distinguishes acceptable designs from 
unacceptable designs, can be grouped into two categories. Any 


* The expressions for the diameter of annulii Ai are 
already obtained in 3.11 and 3.12. 
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constraint which restricts the range of design variables for 
reasons other than the direct consideration of performance of 
the design is called a side constraint. A constraint which is 
derived from explicit consideration of performance of the 
design is called a behaviour constraint. 

SIDE CONSTRAINTS 

Since vector notation is used in the formulation of 
the problem, a set of side constraints must be imposed on the 
magnitudes of the vectors i.e. 

Si ^ 0 i = 1,2, 3, 4 
BEHAVIOUR CONSTRAINT S 

(a) Maximum allowable stress constraint - As discussed 
under 3.25 and 3,26, the dynamic stress induced in the material 
should be less than the allowable stress as 

21 

find = — -2-0 , i = 1,2, 3, 4 

k1T m^ Yi Si 

3.01 

fall = fst 3,01+V 

m-i = B(i) i = 1,2, 3, 4 

m^ = BCl) 

m^ = B(2) 

m3 = B(3) 

ni4 


= B(4) 
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V 


i 


-nDj_N 


FINDj_ = find 


FiJLLj_ = fall 


The constraints are represented as G(I) as given below : 
G(l) = FINDl-FALLl 
G(3) = FINDS -FAILS 

G(5) = FIND4-FilLL4 


As these calculations are for sun gears, and sun gear of 
first and second train is same, the constraint is expressed 
only for the first gear train sun wheel only. 

i'llternativel^.y, using Buckingham equation for dynamic 


load, 


G(l) 

= WBl - 

1.S5 WDl 

G(S) 

= ms - 

1.S5 "WDS 

G(5) 

= V©4 - 

1.S5 \-jD4 


(b) Neighbourhood Constraint 

(Pg). + 2 

G(2) = CIl - 3.0 
G(4) = CIS - 3.0 


1 = 1,3,4 


G(6) = CI4 - 3.0 
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(c) Ann\jlil Size 

The diameter of second, third and reverse train 
annulii is restricted with respect to the diameter of the 
first train annulus as 

D1-D4 

DIFl = i'lBS 

DIF2 = ABS 

G(7) = 0.85 D1 - D3 

GCS) = 0,15 - DIFl 

G(9) = 0,15 - DIF2 


(d) |4odule Co nstraint 

Tha mlninium value of modules used in any of the train 


should he greater than 1 i*e. 
B(i) = m^ ^ 1 


Or 


G(10) = B(l) - 1.0 


G(ll) = B(2) - 1.0 
GCi2) = B(3) - 1.0 


G(i3) = B(4) - 1.0 
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(e) ye],p.clty Eatio Constr^ipt 

The velocity ratio is allowed a variation of 2 per 

cent from the prescribed values. These are expressed as 
follows : 

(VR)^ = B(i) i = 5, 6, 7, 8 
Or 

G(14) = B(5) - 4.195 

GCl5) = 4.36 - B(5) 

G(16) = B(6) - 2.S814 

G(17) = 2.4786 - B(6) 

G(18) = B(7) - 1.5582 

G(19) = 1,6218 - B(7) 

G(20) = B(8) - 5.8306 

. G(2l) = 6.1094 - B(8) 

( f ) Interf eren ce Constraint 

To avoid interference, there is a limitation to the 
number of teeth of sun gears (3.5) i.e. 

Sf, - 16 .^0 

Representing SI = E(9) 

52 = B(10) 

53 ? B(ll) 

54 * B(12) 
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These constrairits are formulated below : 

G(22) = B(9) - l5«0 
G(23) = B(lO)- 16*0 
G(24) = B(ll)- 16.0 
G(25) = BCi 2)- 16,0 
(g) Assembly, constraint 

The assembly constraint is not used in the programme, 
because it cannot be used as such. This is given by 

% + 

= Integer i = 1,2, 3, 4 

All these quantities a 2 ?e integers. 

4.2 SOLUTION SCHEME 

The method chosen here is the unconstrained minimi- 

qo 

zation : Variable. Metric Method of Fletcher and Powell . 

The number of teeth of gears and the number of 
planets must be integers. The modules must be selected from 
the preferred numbers. The assembly constraint consists of 
all integers. These requirements make the problem as Non- 
linear with mixed integer and real variables. 

The only way to solve this problem is to float the 
variables and obtain the minimas. Then these variables are 
given the nearest integer values or preferred values in the 
case of modules. These modified values may or may not lead to 
the minimas. As this method is used only as a check for the 
previously obtained results, the values of variables are compa- 
red with the values obtained from optimum solutions of each set. 
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The scheme of solution of this method is briefly 
discussed below ; 


It requires only first order partial derivatives of 


the function with respect to the design variables. The steps 
required for the minimization of 4 function are as follows : 

(i) To select a starting point and [h°] , where tn^jis 
any positive definite matrix of the order nxnj where n is 
the number of design variables. In general to start with[H<^] 
is taken to be an identity matrix, 

(ii) To find the gradient at x^ , For the i'tli 

iteration x^, the gradient is Vf(x^) and the positive definite 
matrix is |jl\! . 

(iii) Find the next feasible point as follows 

Let = - [h^] X/f (?-) (4.23) 


-i V -i 

and 


x^"^^ = x^ + “Cr 


(4.24) 

(4.25) 

vi 

where a is a scalar and is greater than zero, which is chosen 

"“i “"i 

such that it minimize f along S starting at x . 

(iv) Evaluate y^ = '\jf (x^*^^) - Vf (4.26) 

(v) To solve the following matrix at (i+l)'*"^ iteration : 

( [H^Ky^) (y^)^ [#] 


[h^'*‘1]= [h^] 


( (7^)^ [h^] (7^) 


(4.27) 


(vi) To begin the next iteration from step (iii) 

(vii) To stop algerithm when convergence criterion is satisfied. 
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4,2: •1 CONVERSI ON TO UNCON ST RAIIMBD MTNIMI2.ATI01T 

The general mathematical problem is to determine a 
vector X that minimizes the fimction f(x) subject to the 
constraints 

gi(x) >0 i = 1,2, ,in 

33 

Fiacco and Me Cormick have developed an algorithm for trans- 
forming this mathematical programming problem with constraints 
into a sequence of unconstrained minimization problem'. The 
procedure is based on the minimization of a new function 
Called the Penalty function defined as : 

P (J,r) = f(;) H- 2 7^ 

i=l 

over a strictly monotonic decreasing sequence of r-values, rk. 
The conditions on f(x) and gj_(x) are : 

(a) f(x) and g^^Cx), i = l,,....,m are convex functions and 

(b) P (x,r) is strictly convex in the interior of the 
constraint set for every r ^0»- 

Then, the optimal solution to the constrained problem appi^a- 
ches the global minimum of the constrained problem as r 
approaches zero. 

For r 2 _ 0* and a point x ° strictly within the 

constrained set that is gj^Cx '^) i = Ij ••*•••■ jia, it is 

expected that a miniiaum of P (Xjt) exists inside the cons- 
trained set, since on the boundary of the constrained set 
some gij^(x) must be equal to zero and P^Cx^r^) depends on the 
valne of r^j and is denoted by xCr^)» The point xCr^^) lios 



48 


Value of r^, and is denoted byx(r]_). The point x(r^) lies 
inside the constrained set« By reducing the value of nj the 
influence of sumniation temij which penalizes the closeness to 
the constrained boundaries, is reduced, and in the minimi- 
zing P^(x,r), more emphasis is placed upon decreasing f (x) • 

If the minimization process is repeated for decreasing sequence 
of r values, each minimum point xCr^j.) can move closer to' the 
boimdaries of the constraint set if it is profitable in 
terms of reducing f(x) . This method is particularly attrac- 
tive for the problems which have non-linear constraints, 
since it approaches the solution from points inside the 
constraint set and thus avoids the difficult movement along 
the non-linear boundary of the set. 

If 

(a) the interior of the constraint set is nonempty 

(b) the function f(x) and -g^^Cx) , i = 1,2, ,m, 

a". 3 continuously differentiable 

(c) the set of points in the constraint set for which 

t(x) is bounded for every finite Vq, and 

(d) the function f(x)is bounded below for x in the 
constraint set, then 

there is optimal solution to the constrained problem as, the 
value of r approaches zero. 

Now to apply the Fiacco and Me Cormick algorithm to 
this problem, the Penalty function is formed for the gear box 



'4a 

volxime :• 

P^(x,r) = V + rZ — T-T- 

" Ci(x)' 

r 1 + , 1 ^ 1 

^ ^ 1 FIMDl-FALLl Cli-S.U Fin)2-FALL2 

1 . 

CI3-3.0 "* 

The conditions mentioned above in this section, if they are 
not satisfied, it is possible that the algorithm may not lead 
to a global minimum. In such a case one has to start from 
different feasible points and obtain the minimum of the 
minima* If on the other hand minimum obtained from different 
starting points is same, one can conclude that the minimum 
obtained is also global. 

4 ,2 *2 Computati onal A snect of Penai.tv Function 

Following steps are involved in this technique : 

(i) To choose an initial feasible point Xq . 

(ii) To select r^, the initial, value of r. The selection of 
initial value of r is quite important. Initially it 

would seem advantageous to select r., = 0. where ^ is 

selected as small as possible, without loosing significance 
due to round off errors, etc,, but it effects the rate of 
convergence . 

Below is given a method of choosing r^, which is 


usually adopted 
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Since a necessary condition for (I,r) to bo a miniamm 
is the vanishing of the first partial derivative a natural 
choice of r-L would be given by r that minimizes the magnitude 
of the square of the gradient of P^Cx,r) at i.o. 


2 

P^(x,r) j = min^ [ y f(x) 


V gi(^)!g2(^) 


(4*29) 


For simplicity, let P(x) = n' g , then r^ is given by, 


P I pCf) 


— o V 2 


(4 *30) 


Equation (4,29) is obtained by differentiating (4,28) with 
respect to r and equating the resulting expression to zero. 
Equation (4,30) can be used if r^N 0 which is obtained if 

V p 40 . If V' p(^) >0 (4.31) 

gives ^ decreased at its 

maximum spatial rate, without increasing the the 

P-function, by proceeding along - V f (^) . In this case, 
therefore, taking a sequence of steps of given length down 
the gradient f and recomputing (4.30) at each new point until 
a positive r^ is obtained, is an efficient way to proceed, 

(iii) To determine minimum of Pf(x being the value 

of r at the Iteration. 
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(iv) To terminate solution of the above step if the conver- 

Lgence criterion is satisfied* u. convergence criterion is 
when j V'' P^(x)j , where 0 , but a small specified 

quantity* 

(v) To select = *^i ^ 15 ; where 0 'C 5^^,^ and to continue 

the procedure from step 3 for the new value of r. 

(vi) To stop algorithm, when the final convergcnco criteron 

is satisfied, i. final convergence criterion is one when 
m 

r ^ 1— - where S \0 is a very small quantity. 

gi(x) 

For the Fletcher & Powels method, the gradient of the 
objective function and constraints at x are required. These 
are obtained by finite differences method as ; 


4.2,3 .The Gradient 

The gradient of a n-dimonsional function is defined as : 

(4.32) 


T!e) = vFa) = (^,^, Hr) 


n 


Or it may bo expressed as 

m 

V 0(x) = r (4^ + r 2 


Ogi 




+ r 


'> 1 >gi 


i=l gj g'i 

m m 

* 'S -i iflU 

* • * • , C ” r p M 


■bx. 




(4*33) 


tlT. KANPUR, 




of the interior penalty function for 


Thus the gradient 
a fixed value of r is 


0 U) = (-^ , 

liXi -^Xg 

or it may be expressed as 

\ 70 (x) = [(in + r '> V . -^) (-^ 
I 5^1 T=i ef 


(4.33) 


m 

'=^7 1 gi 

i=l ^ ^ ^2 


( 




m 


■bx, 


n 




1=1 g 




(4.34) 


In most design problems of even moderate complexityj it 
is simply impractical to derive analytical expressions for all of 
■the £ and — ^ terms. Therefore it is necessary to resort to 

JfL 

finite difference derivations. 

The difference between the accuracy of eq,(4,33) and 


that of eq.(4,34) is due to the differences among the functions 
0, f and g. When x is near a constraint surface and r is small, 
0 becomes a very distorted function while f and g usually remain 
well behaved. Examining the Taylor series expansion of a 
function of one variable about the point x^. 
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Solving for the first derivative yields 

_ / \ 1—1 / -N 1 ^ 1? \ / / 1 






(4.36) 


If the second and higher order derivatives can be 


neglected, then 


ax ' 




(4.37) 


is a good approximation to eq,(4.36). 

The gradient obtained from eq.(4.37) and (4.34) is 
accurate enough for use in the variable metric method. 

4.2.4. Determination of sten Length 


-1 

X is given by ; 


If x^ is the initial feasible point, the move to reach 


. ^[5^1 (J°, r) 


( 4 « 38 ) 


where (^,r) is the gradient vector and a, 0 is the 

step length chosen to minimize along - [h^J x^P^C?^), so that 

— 1 ^P 

X is a feasible point. ,, is determined as follows : 

\p 

First an upper limit 7 ~. is determined so that none 
of the non-negativity restrictions are violated. Let the 
vector - [h^] ^ V Ff(^j i") is denoted by 


=[a^i = .[ll]?=’VPf(?’,r) 


let fi 


- ^ i > i - 


(4.39) 


(4.40) 



^ i.Pthen none of the non-negativity restrictions are 
violated. 

First this ^ is divided into t>ro equal parts and more 
gradually do-wn each time tahing a step length equal to 0,5^ , 

At each new point- it is checked that along this direction, 
whether the P-function has decreased and. the new point is the 
feasible point. 

If all the constraiJits are satisfied and the P-function 
has decreased from its previous value in that direction, then we 
further move in the same direction. If any of the above 
conditions are not satisfied, the step size is reduced further 
to half. If the step size becomes a small fraction of p , 
consistent with the required precision, say one thousandth of ^ , 
the move along that direction is stopped. Thus is deter- 
mined by this numerical search. 

After getting the point from x°, further x^, x^,.,*. 
are obtained by application of the same procedure. The P- 
function is approximately minimized for a particular value 
of r at a point when P^Cx^"^^) - Pj.(x^) = 0 

The r value is decreased by quarter and half and the minimization 
is continued for a decreasing sequence of r. The algorithm is 
stopped when the change in x minimizing the P- function for two 
successive values of r is less than 0#0000l and the difference 
in the value of objective function and P-function is also less 
than 0.0005. 
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i 



Return 






5S 
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SUB RODTLUS GRAD 


Set x=0.0005! 


] 


pall OPUi'lC 


Y 


i Set Z2=W 
X3=P 


1 


Confute 


'HDUte V I 

b'(i)=b(i)+ a i 

t 


Call OFDl'IC! 

f * 



Repeat 11 times 
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The results obtained from penalty function approach are 
listed helov* 3ach set of variables found^ correspond to a 
feasible solution. 









CHATTER V 


DI^US SIOH of rb stit.th 


The results of the computer aided design are given 
in Table I of Chapter III, Ey taking different numher of teeth 
on the sun gear of first gear train, the second, third and 
reverse gear sizes are obtained by using iterative procedure. 
From the various .solutions obtained corresponding to each of 
16-36 number of teeth of sun gear of first gear train, the 
minimum volume solution has been picked up from each of the 
subsets. From these minimum, solutions, the minimum most has 
been selected as the best design. 

While formulating the design procedure of gears, 

Barth equation and Buckingham's equation were used separately. 
The Barth equation, being less conservative, leads to much 
lower values of volume as compared to the Buckingham’s equa- 
tion. The lovj-est being approximately 52?^ of the existing one. 
The number of solutions also were much more than those obtained 
from Buckingham's equation for the same tolerance on the 
diameter of annuli of the gear trains. In the present case, 
the design has been based upon the Buckingham’s equation, 
because it predicts nicely the value of the dynamic load and 
it has given more satisfactory designs in practice. 
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In order to check that the optimum solution by the 
search technique is actually the minimum most, the penalty 
function approach is tried. Some of the minimum solutions 
of the search technique were fed as initial feasible solutions. 
After somo iterations it leads to almost the samo solution, 
which shows that the initial solution was near to the local 
minima. The results in different iterations also show that the 
slight variation in the variables violates one or the other 
Constraint indicating that the feasible region is very small 
and the penalty function is giving the local minima. In other 
words, the possible design space, in the case of epicyclic 
gear drives is not a continuous one. Bounded by constraint 
surfaces, it constitutes rather a number of discrete patches 
in the design surface (mainly because of tight constraints on 
the number of teeth). 

T-he optimum solution (Design No. 7) obtained occupies 
about 78 percent of tho existing volume of the gear box which 
shows that there will be proportionate saving of material, 
space occupied and all other dependent factors. From the 
various solutions obtained, one of the solutions (Design No. 16) 
is almost exactly the same as tho existing one which justifies 
the selection of material, strength and the procedure adopted. 



APFEKDIX A 


^ SCRIPTIOH OF SOME CRIT ICiJ. PJJffS /ilTO PRD-RT.TiM.q 


The gear box provides a means of overcoming the 
initial inertia of the automobile and assists the engine "while 
it is in motion by making the minimum demand on its powerj 
when extra resistance is encountered. The gear box also provides 
a means of uncoupling the engine from the remainder of the 
transmission when used in conjunction with the clutch, 

M/iTSRI/iLS 

The pressure on the teeth of gears is extremely high 
and the material employed in the manufacture of the gear must 
possess toughness and hard wearing qualities. Usually a high 
nickel content is called for and such materials as nickel- 
chrome molybdenum steel and nickel molybdenum steel are used. 
They are heat treated in the course of manufacture of the gear. 
The same mr.terial is also used for the shafts. This also faci- 
litates pi lions to be integral with the shaft. In the forma- 
tion of cauing and extensions, cast iron is commonly used, but 
some manufacturers favour aluminium alloys on the score of 
lightness and die castings are fairly common. Bushes are of 
phosphor bronze or lead bronze. 

The Wilson Gear Box foims the basis of all the 
different makes, as they only vary in detail such as selector 
and toggle action, oil pump or lubrication and design of 
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clutch for top speed. All types retain its outstanding features 
i.e* employ servohralces in such a manner as to cancel out the 
braking strains on the annuli and the compounding of different 

trains so as to reduce speed and stresses of 
planet gears. 

DESCRIPIIOH OF SOl'E PARTS 
Brakes and clutches 

The important advantage which epicyclic gears have 
over layshaft gears is the ability to change from one ratio 
to another, without loss of torque transmission. This is 
achieved by suitable clutches or brakes associated with the 
reaction members, one of which is released as another is 
engaged. This is called ‘Hot Shift’. 

Since a torque balance must exist when transmitting 
power, it will be apparent that the torque in the brake or 
clutch holding a reaction member is the difference between 
input and output torque. Thus for large gear steps, the reac- 
tion torque is relatively high and it is common practice to 
use a self-energizing brake. For lower ratios where the 
reaction torque is small, plate or cone friction clutches are 
more commonly employed. The potential speed of a reaction 
member is inversely proportional to the torque required to 
hold it stationary i.e. a reaction member, which by virtue 
of the gearing, would run at high speeds with the output 
member stalled, requires only a small torque to hold it 
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stationary* The actual energy which is absorbed (and dissi- 
pated as heat) in a friction clutch or brake during a shift 
is related to the torque transmitted and the ratio step and 
is independent of the potential speed difference. 

Bands are useful for holding high tcrque reaction 
members as, for Instance, with low forward and reverse drive# 
The self -wrapping servo effect of the band can be used to 
advantage in such cases, but band brakes of this type may be 
difficult to control to give smooth engagements and may not 
readily release on disengagement. 

Band brakes are usually applied by hydraulic servo- 
cylinders acting directly, or through levers, on one end of 
the band, while the other end is anchored to the casing. 

Care should be taken to avoid high loadings being applied to 
the casing at the servo-cylinder and the band anchorages# 

Plate clutches used are usually of the multi-d^sc 
tvpe » The^o give a large area of surface and, hence, energy 
capacity. These clutches are engaged by axial forces which 
are derived from springs or hydraulic cylinders. For dynamic 
engagements, smoothness is improved by modulating the hydraulic 
pressure. To ensure quick disengagements and so avoid torque 
reversals, it is usual to separate the clutch plates by 
light springs f 

Cone clutches are not widely used with epicyclic 
gear trains because the energy which a cone clutch can 
absorb during an engagement is limited as compared with the 
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multidisc clutch, hut other important advantages are available 
smaller space, particularly in axial length, lower cost, lower 
clamping load due to wedge section of the cone and lower drag 
when disengaged,. 

Bearings 

In considering the bearings for the elements of an 
epicyclic train, the requirement for supporting the main 
elements is simplified by the inherent balance obtained by 
the use of two or more planets. Concentricity is important 
to ensure equal load sharing at each tooth contact point. In 
the absence of perfect concentricity it has been found that 
considerable freedom of one or two of the elements can permit 
automatic centering under load. 

The most, difficult bearing problem in an epicyclic 
gear train is that of the planet pinion bearings. The 
problem can be more apparent in an overdrive unit. At high 
road speed, the transmission is normally 1^1 ratio with gear 
trains rotating as a whole, but with an overdrive the planets 
are working continuously over very long periods. The planets, 
as well as turning on their own axes, orbit about the axis of 
the carrier, thus developing very appreciable centrifugal 
force at high carrier speeds. For this reason planet bores 
must be as large as possible to reduce weight, although from 
a consideration of tooth bending stresses, the amount of metal 
below the tooth root should not be less than the full tooth 
depth. 
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Th0 UGodlG boaring has come to he widely used for 
the planets of opicyclic trains and although some lubrication 
is essential, the necossity for full pressure oil supply does 
not arise as with plain bearings. ■ 

Making the gear bore as large as possible permits the 
use of a pin of minimum deflection. This is important with 
needle bearings since any appreciable deflection causes locally 
high loading on rollers.. The distribution of loading on the 
bearings, even of simple planets, is not uniforn owing to the 
effect of thrust couples and centrifugal force, it is therefore, 
usual where space permits, to use two separate sets of rollers 
suitably spaced apart. 

Moreover, the centrifugal force acting on the needles 
becomes a problem as a result of the tooth loading and the 
centrifugal force moving towards the underside of planets. The 
needles which are loaded on the underside roll as they should, 
but where tie maximum clearance exists at the point opposite 
the resultant load, the needles tend to accumulate and lose 
their rotational speed. When passing into the load area 
again they skid as they are accelerated to full speed. 

If there is excessive diametral clearance, loose 
rollers can run out of parallel with the planet pin axis, 
causing further skidding and heat generation. For this reason, 
caged roller bearings are sometimes used, where a smaller 
number of rollers is employed, with a light steel or aluminum 
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cagG to hold tneni parallel with the piris. Such caged tGarings 
themselves introduce fresh problems, by virtue of the centri- 
fugal loading on the cage itself. 

All these problems are associated with high speed 
rather than with high load. Appreciable advantages can be 
obtained by improving the lubrication over that which obtain 
with simple oil grooves, and one successful method is to centri- 
fuge oil from the hollow main shaft, collecting it in a labyrinth 
arrangemont and guiding it into the hollow planet pin and out 
of a suitably placed hole into the roller bearing, 

CASING 

The design of the casings for epicyclic gear boxes 
is relatively simple since tooth loads are balanced and bearing 
loads on the main elements are nominal* For this reason it is 
possible and usual to use aluminium alloy casings and, on 
transmissions intended for large volume production, these 
should invariably be designed as pressure die castings* 

NOISE 

The principal source of noise is deformation as teeth 
enter and leave engagement, A study of gear noise reduction 
should concentrate on gear design to reduce noise at source, 
but it should be remembered that resonance of some other 
part of the vehicle may magnify, to unacceptable levels, 
noises which are small in themselves* 
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almost always difficult to diagnose and cure. 
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SPBEDS IK Ml EPTCYCL IC TRAIH 

The various possible coEibirLatious in an epicyclic 
gear train are described belou with the help of an example. 

1, In most Cases, the sun gear A is driven from the engine and 
the carrier of the planet gear D is the output member conveying 
the drive to the next transmission member* The outer cylindri- 
cal surfac' of the annulus gear is provided -with a band brake, 
such that when this bralce is applied, the annulus gear is 
locked. Assuming that the annulus (B) is locked and power 
applied to the sun gear Fig. Ci), then the carrier (D) will 
revolve, carrying with it the planet gear (C), If the gear 

(A) rotates in a clockwise direction, (C) will rotate anti- 
clockwise and the carrier (D) will move in a clocki^ise direc- 
tion also. Therefore if the carrier is a power output member 
it will rotate in the same direction as the sun gear but at 
a slower speed than the sun gear. 

The speed of rotation of the planet carrier can be 
A 

shown to be times the speed of the sun gear (A) where A 

and B denote the number of teeth of gears (A) and (B) . 

2, If the carrier CD) is locked, as by a clutch unit, and the 
sun gear is power-driven member, then as shown in Fig.(ii), 

if (a) revolves clockwise, it will rotate the planet gear (C) 
anticlockwise and this gear will in turn rotate the annulus 
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gear anticlockwise and at ^ times the speed of U) . 

3. If the sun gear (A) is locked and the carrier (D) is driven 
in a clockwise direction Fig.(iii) then the annulus gear (B) 

driven hy the planet gear (C) in a clockwise direction 
and at a higher speed than (A), the ratio of the annulus gear 
to carrier speeds being 

4, If the annulus gear (B) is locked and the carrier (D) 
becomes the clockwise direction drive Fig..(lv) it will rotate 
the planet gear C in an anticlockwise direction and s-un gear 
(a) in a clockwise direction and at a higher speed than its 
own. Thus the sun gear will rotate at times the carrier 
speed * 

5. If the Carrier D is locked, the annulus be po^-:er driven 
and the drive is obtained from the sun gear which will rotate 
at ^ times the speed of (A). 

6, Again 3.f the sun gear A is locked and annulus gear B is 
the driving gear, the carrier D will rotate at a slo'v/er speed 
given by times that of the speed of B. 
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1 .'^F PLAMLTS:!-,5X,- 
PP IMT20 
F1=6300o0 


EETH^sTX 9 -^ANNULUS TEETH-^f OX? ^PLANET 
FAp PAT 10* » 6X )«-MODULE-'<- s 6X 5 *ANMULUS 


TEETH*>3X » *N0 
DIAMETER* ) 


RPM=1100o0 

Ti= 62 o 1 , ' 

bK= 3.5 

F.ST = Fl*l00oy/7o0 
AC=0o003? 

AX = “0.. OOOO/lOCOoO 
P I =Ao*ATAN( 1 o ) 

Xr 1 = Ac ?8 

THE NUMBER OF TEETH ON SUN GEAR ARE MADE TO VARY FROMlS UPtfAKDS 
^0 A5 I = 1 9 20 

^ 1 = 16 qO+AI ' 


calculation for the number OF TEETH FOR ANNULUS ANCPLAI-TUT GEARS 
Al=(XRi-l,0)*Sl 
Pl=(XRl-2oO)*Sl/2oO 

THE number of TEETH ARE ROUNDED TO THE NEAREs T I NTEGEr VALUES 
NS1=S1 


NAlxAl+0 = S 
NPl=Pl+0.5 
SN1=NS]. 

AN1=NA1 

PN1=NP1 

CALCULATION FOR LEWIS FORM FACTOR 
Y] =0c 15A-0„012/SN1 

NEIGHBOURHOOD CONDITION. APPLIED TO CALCULATE MAXIMUM NUMBER 
Kl=(PNl + 2oO) /(SNl+PNl) OF PLAl-ISTS 

THETA1=AR5IN(R1 ) 

PMAX1=PI /THETAl 
NPMAX1=PMAX1+0 c 5 
PMAXA1=NPMAX1 
Do 40 N^l>6 
BN=N-1 

PNX(N):^PMAXA1-DN 
IF (PNX(N ) oLToPoO) GO TO AG 

check for ASSEMLlY OF GEARS IN THE FIRST GEAR TRAIN 
FCT0R1= ( ANI+SnI ) /PNX ( N) 

1FAC=FCT0R1 

fac=ifac 

IF {FCTOR1-FAC-0o0000001)259A0»40 
NP-PNX(N) 

PTl=NP 

selection of module FROM THE PREFERRED VALUES TO SATISFY 
I'-iAXliiUivI ALLO'IA'iLS STRESS 
Do 90 M=l>10 
A|M = m 

Am 1 = AM/2 oO + l oO 
DIAS1=AM1*SN1 
Dl AA1=AM1*AN1 
IF(UIAA1 oGTc260o0)G0 TO 90 


C 

C 
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70 

75 

77 


VEL ‘il-nl '^DIA.0^'M^p^<,/9s.4/l?o0 

vn 'fl 1/01 ASi^f 1 0U0,0-x-2c :>/PTl 
B1 = PI#AN'1'«Sk/25o4 ' 

Fp'-] VnSl+AC 

CU-166O.U*Lrpi/0.O0i 

WH 1 = F ' T if B 1 'i Y 1#P I /2 ,4 

format (10X,4(UX, 15) ,3F16o?) 

PR I NT 77 
Format (.? ( i:-;- ) ) 


calculatiomo for the size of second gear train 

Do Z50 L=l»10 
al = l 

XRZ-^ZoOP + Al/IOOc 

BN2=SNT 

NAZISM? 

CAL aJ|./\T ION FOR THE. NUMBER OF TEETH FOR ANNULUS ANDPLAJTET G3ABS 
A2= ( .sNl+‘XRT,,'t ( Io0-XR2 ) ) /(XR2-XR1 ) 

PT-" ( A2-.SN2 ) /Z„0 

THE NUMUER OF teeth ARE ROUNDED TO THE NEARESTINTEGER VALUES 
NA2=A2+Oo5 
NP2^P2+0, 5 
A,N2 = NA2 
’30 PN?=NP? 

calculation for lewis form FACTOR 
Y2 = 0. ] A4*-0«qi2/SN? 

NE IGHnOURMOOU CONDITION APPLIED To CALCULATE MAXIMUM NUMBER 
RP--^ {PN2 + 2eO) / (SM2+PN2 ) OF PLANETS 

TIIETAZ^ARSIN (R? ) 

PMAX2 = RI /T1IETA2 
NPMAX2';=rMAX2 
r’MAX2 = NP;iAX2 
DO 2AU NI'-UO 
IjN'-NR-i 

P2X (N ) :-.PMAXr-ISN 
IF(R2X(N) oLT.. 2.!;;)GO rc 24C 
CI-iFlk FOR' ASCLLDI.Y OF GEARS IN THE FIRST GEAR TRAIN 

Ft,. mloi'+'.i.iz ) /p:ix (.N) 

IFAC-FrTOR2 

FA(=:IFAC 

IF { FCTOP2-FAC-0.0000U01 ) 225,240 , 240 
25 PT=P2X(N) 

NPT=PT 

PT2=NPT 

selection of module from the PREFERRED VALUES TO SATISFY 
hAXIhm AU .0 I/iPLE STRESS 
Do 290 MAbUIO 
tJM = MA 

AM2=Bm/2oO+ 1 .0 
Dl AS2 = AM2^fSM2 
DIAA2=AM2«AN2 

THE Diameter of the second annulus is taken nearly equal to 
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c THE Diameter of the first annulus for cc-'pactness 

IF (DIAA2oGTo0c>99*DImA1 oAnDoDI AA2cLTc l = 01-»-DI AA-- )Go TO 235 

Go TO 290 

0 check For strength of gear tooth due to DYNA^IIC load using 
c BUCKINGHAM EQUATION 

2 35 VE[_S2 = pi*DIAS2*RW/25o 4/ 12cO 
■' W2 = 2»0*T1/DIAS2*RPM^^2o2/PT2 
B2 = PI*AM2*SK/25 oA 
ERR2=AX*VELS2+AC 
WB2 = FST*B2*Y2^f-P I*AM2/2 5oA 

WD2 = 0 <.05->f\/ELS2* ( B2*C1 + W2 ) / ( Oo 05*VEL5'2-rSORT ( B2*Cl-i-W2 ) )+W2 
I F {WB2-1 o35-«-WD2 ) 2 90 , 2 70» 270 
270 PRINT2755NS2jNA2 ,MP2 , NPT , XR2 > AM2 , D I AA2 
275 Format (10X54(11X,I5),3F16o2) 

290 Continue 
2A0 Continue 
250 Continue 
PRINT277 

277 Format { 2 ( iHO ) ) 

Do 350 MC= 1,6 
XR3=1o56+AM/100o 

c ^calculation of the number of teeth of annulus and planet gears 

G11=AN2*( ANl+SNl ) -AN2*SN 1*XR3 

G12= ( ANl + SNl ) *( AN2+SN2 ) -XR3-x-ANl*SNl-XR3-:HSL'l«-A:-i2-XR3*SN'' *SM‘' 

G3=G11/G12 

Do 345 NC=1,10 

C THE number of TEETH ON SUN GEAR ARE MADE TO VARY FR0M16UP./APDS 

AN=NC-1 
AM = MC 

B3=16o0+An 

C ^CALCULATION OF TH5 NUMBER OF TEET8 OF ANNULUS AND PLANET GEARS 

A3=S3/ (G3-lo0 ) 

P3=(A3-S3)/2oO . 

C THE number of TEETH ARE ROUNDED TO THE NEARES T I N TEGER VALUES 

NS3=S3 
NA3 = A3+0<.5 
NP3=P3+0o5 
SN3=NS3 
PN3=NP3 

IF(PN3oLEc,15oO)GO TO 345 

SN3=N53 

AN3=NA3 

C calculation FOR LEWIS FqR*'' 'ACTOR 

Y3=Oo154-Oo912/SM3 

C neighbourhood CONDITION APPLIED TO CALCULATE MAXIMUM NUMBER 

R3=( PN3+2oO) / (SN 3+PN3 ) OF PLiNSTS 

THETA3=ARSIM (P3 ) 

PMAX3=PI /THETA3 
NPMAX3=PMAX3 
PMAX3-NPMAX3 
Do 340 KA=1,6. 

Bk = KA 

PX (KA)=PMAX3-6K 
IF(PX(KA)oLTc2o0)G0 to 340 

C check for assembly of gears in the first gear TRAU! 
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FCT0R3= ( AN3+SN3 ) /PX (KA) 

IFAC=FCT0R3 

fac=ifac- 

IF (FCT0P3 -FAC-0 cOCO 0001 ) 3259340,340 
325 NT3=PX(KA) 

PT3=NT3 

selection of module from I>E preferred values to satisfy MAXiMUf 
Allowable stress 
Do 390 MD= 1,10 
Dm=MD 

AM3=DM/2oO+ 1 oO 
DIAS3=AM3*SN3 
Dl AA3=AM3*AN3 

IE(DlAA3oGToOo90^DI/-\AloANDoDIAA3oLl'^jOop3'if'DIAAilGo TO 335 
Go TO 390 

THE diameter OF THE FIRST ANNULUS FOR COMPACTNESS 
BUCKINGHAM EQUATION 
3 35 VElS 3 = PI*DI AS3*RPM/2 5o4/ 12o0 

W3=2o0*Tl/DIAS3ia000o0'2»2/PT3 

B3=PI*AM3*SK/25o4 

ERR3=AX*VELS3+AC 

WD3 = 0o05*VELS3* ( B3*C1+W3 ) / ( 0 a 05"-f-VELS3 + 5QRT ( B3i:-Cl+W3 ) )+W3 
WB3 = FST*B3*Y3^tP I*AM3/25o4 
IF (WB3-1 o35^<-WD3) 3 90 9 370,370 
370 PRINT375>NS39NA39NP39MT3 ,XP 3 » AM3 , l I AA3 
375 ForMAT(10X94(11X,I5) ,3F15o2) 

390 Continue 
340 Continue 
345 Continue 
350 Continue 

PRINT377 

377 FormAT(2(1HO) ) 

C calculations for the SIZE OF REVERSE GEAR TRAIN 

Do+145 ND= 1,20 

DN=ND-1 

Sr=20oO+DN 

c ^calculation of the number of teets of annulus and planet gears 

AR = ABS(SNl*SRrr(XR4-loO)/ANl ) 

PR= ( AR-SR5 /2oO 

C THE number of TEETH ArE RC-UNDED TO THE NEAREST INTEGER VALUES 

NSR=SR 
NAr = AR+0 o 5 
NPR=PR+0o5 
Snr=nsr 

ANR^NAP 

PNR=NPR 

IF(FNRoLEo15o0)GC TO 145 

c calculation for lewis Form factor 

Yr=0o 154-Oo912/SMR 

C neighbourhood CONDITION APPLIED TO CAlCULA TB 4A}CIMUi4 IiUl-'IBER OF 

RR=(PNR+2o 0) / (SNR+PNP.) PLjU^IEIS 
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THETAR=ARSIN(RR) 

P,'''*AXR = PI /THETAR 
NPMAXR=PMAXR 
PMAXR=NPMAXR 
Do 140 KC=l,6 
CK=IAC 

PRX(K)=PMAXR-CK 
IF(PRX(K) oLTo2oO)GO TO 140 
FCT0RR= ( ANR+SNR ) /PRX ( 2 ) 

ifac=fctorp 

FAC=IFaC 

IF (FCToPR^FAC-OoOOOOQOI) 125,140,140 
125 NTR=PRX(K! 

PTR=NTR 

c selection of module from the preferred values to satisfy maxim 
c allowable stress 

Do 190 MR= 1,10 
RM=MR 

AMR=RM/2oO+1CO 
Dl ASR=AMR*SNR 
Dl AAR=AMR*ANR 

IF {DIAaRo6To0o99^"DIAA1..anDoDIAARoLTo1o 01*DI AAi )G0 TO 135 
Go T0190 

C check for strength of gear tooth due to dynamic LOAD USING 

C BUCKINGHAM EQUATION 

135 VELSR=pI*DIASR*RPM/25o4/12o0 

WR=2oO*T 1/DI ASR*1000o0*2o 2/PTR 

Sr=pi*aMR*SK/25o4 

frrr=ax*velsr+ac 

Cr=1660o0^ERRR/0o001 

WDR = 0o 05*VELSR* (BR*CR+WR ) / ( 0o05i: VEL5P+SCRT (BR*CR+WR) ) +WR 
WBR=FST*BR*YR*P I*AMR/25o4 
IF (WBR-1 c 35*WDR ) 190,170,170 
170 PRINT175',NSRTNAR,NPR,NTR , XR4 , AMR , D I AAR 
175 FORMAT! 10X,4( IIX , 15 ) ,3Fl6o2 ) 

190 Continue 
140 Continue 
145 Continue 
150 Continue 

print 177 

177 Format 12 ( iho ) i- 
90 Continue 
40 Continue 
45 Continue 
STOP 


SENTRY 
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^ Programme for penalty function appr^'ac-' 

Dimension b(i2) 

Common/ APT/ PI sT , te>fs ,xri ,xr2 »xro »xr4»fend»an 
PI=A oO>ATAN( loO) 

FEND = 63o 0 
AN=1100o0 
T=62100o0 
TE=1o5*T 

read ioo» (B( n , i = is 12 ) 

100 FORMAT! l2F6o2) 

CALL FmIN(B) 

STOP 

End 

SIBFTC FMIN 

Subroutine fmin(B) 

Dimension b ( 12 ) jG (25 ) ^gw r is ) sGp ( 25 ) »delx ( 25 ) sGWi ( 25 ) »gpi ( 25 > >s( : 
1»Y(25),HY(25),H(25»25)5YK(25,25I,AE(25,25)9Y+(25,25) 9HYH{25,25) 
2U2»12)»GF(12>12)»GF1(12) ,DEL(12) jDELKIZ) 

COMMON/APT/PI 5T,TE9FSjXR 1,XR2 »XR3 » XRA > FEND , AN 
COMMON/BPT/ANl 9 AN2 9AN3 ,AN4»2P1 9 ZP 2 9 ZP 3 ,2P4 
7 FormAT(1X»13£9o2) 

80 ForMATi* penalty*, E lO. 2 ) 

108 FormAT( 20X,*THE PRESENT VALUE OF R=*,El2o6) 

102 F0RMAT(1X,12E10o4,2X) 

101 ForMAT(20X,*THE INIATIAL value of R=*,E12o5) 

100 FORMAT! IX, 10 ( E10o4,2X ) ) 

107 FORMAT! 20X,*ThE NUMBER CF ITERAT I0NS=*, 14 ) 

iiu Format (/ lox ,*THESE are the values of constraints*) 

112 ForMAT(/10X,*THES are the values of VRIASLES object function AN! 
IEnALITY function*) 

104 F0RMAT(/10X,*TmESE ARE THE VALUES OF GRADIENTS*) 

109 ForMAT!/10X,*THESE ARE THE VALUES OF DIRECTIONS*) 

38 Call ofunc! b,w,g,p) 

10 Format (//iox,*ani=*,f3«3) 

11 F0R.MAT(//10X,*AN2=*,FSo3 ) 

12 Format ! //iox,*an3=*,F8c3 ) 

13 F0RMAT(//10X,*ANA = *,FF.o3) 

14 FORMAT(//lOX,*ZPl = *>.Fpo3) 

15 FORMAT! //10X,*ZP3=*,F8o3 ) 

16 F0RMAT!//10X,*ZP4=*,F8o3) 

PRINTllO 

PRINT7,(G(I ) ,1=1,25) 

PRINT 80, P 

printio,ami 

PRINTll, AN2 

PRINT12»AM3 

PRINT13,AN4 

PRINT149ZP1 

PRINT15,ZP3 

PRINT16»ZP4 

IF (PeGTo 10o**24 ) STOP 

Call grad(B,gw,gp,G) 

SUMl=0o0 
SUM2=0 o0 
DO 31 -1 = 1,12 
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^UMl-SUMl + GP ( I )*GW( I ) 

:5 i '^UM2-5|JM2 + GP ( I ) -K-GP { I ) 

R=-SUMi /SUM 2 
lF(RoLT=0a0)R=0o5 
PRINTIOSsR 
IR(R) 32*33,34 
32 R=Oo 

F=W+r*p 
Fi = f 

PRINT102, (GW( I ) ,1=1,12) 

Do 35 1=1,12 

35 DELX( I )=-GW( I ) 

PRINT102, (DELX( I ) ,1 = 1 , 12 ) 

SUM=OoO 

Do 36 1=1,12 

36 Oum = SUm+DELX ( I ) *-‘;-2 
Do 37 1=1,12 

37 DElX( I )=DElX(I )/5qrT(SUM) 

OAlL STEF(E,F,G,Fi,DELX,ALP,R,W> 

PRINTi03,ALP 

Oo To 38 

3 3 PRINTlOO, (B( I ), 1 = 1,12 ) ,W,F 
Oo To 50 
34 PPINT101,P 
KE=OoO 

23 Do 39 1 = 1,12 
F=W+r#p 

39 DEl ( I ) =GW ( I ) +R*GP ( I ) 

Do 41 1=1,12 
Do 41 j=l,12 

41 H(l,j)=0c0 
Do 42 1=1,12 

42 . H( I ,i ) = l,o 

N = 0o 

22 YHY=OoO 
STY=OoO 
Do 43 1=1,12 
Do 43 J=l,12 
HY( I )=OoO 
YH(I,J)=OoO 

43 HyH(I,J)=OoO 
Do 44 1 = 1,12 

44 DELX(I)=OoO 
DO 45 1=1,12 
Do 45 J= 1,12 

45 DElX( I )=DElX{ I )-H( I ,J)4DEL( J) 
PRINT112 

PRINTlOO, (B( I ) ,I = 1,12 ),,w,F 
PRINT104 

PR.INT102 , {DELX( I) ,1 = 1,12) 

PRINT109 

PRINT102 9 (DEL(I ) ,1 = 1,12) 

SUM=OoO 
Do 49 1=1,12 
49 Sum=SUM+DELX ( I )**2 
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Do 51 1=1,12 

51 ‘^E3X( I )=DELXUI ) /509T (SUM) 

F1 = F 
F2 = F 

Call step (b,f,g ,fi,delx, alp, r ,w) 

PRINT103 »ALP 

103 ForMAK lUXs-K-THE STEP LENGTH=-;t , E 1 2 o 6 ) 

Do 56 1=1,12 

56 S( I )=AlP*DE1X{ I ) 

Call grad(b,gwi,gpi,g) 

Do 57 1 = 1,12 

DELK I )=GW1( I ) + R^GPl{ I ) 

57 Y ( I )=DEL1 ( I ) -DEL{ I) 

Do 58 1=1,12 

5 8 STY = STY+S( I )MY ( I ) 

IF ( ALP-OoCOOl )67,67 ,51 

61 IF (SQRTtSUM) »LE o 0„01 )GO TO 67 
Do 59 1=1,12 

Do 59 J=l,12 

59 HY( I )=HY( I )+H( I ,J)^f-Y( J) 

Do 60 1=1,12 

60 YHY = YHy+Y ( I ) -ii-HY ( I ) 

Do 62 1=1,12 

DO 62 J=l,12 
YK( I ,J)=Y( I )*Y{ J) 

62 AE{ I , J)=S( I }MS( J ) /STY 

Do 63 1 = 1,12 

DO 63 J=l,12 
Do 63 K=1,12 

63 YH ( I , J )=YH( r , J)+YK ( I ,K ) *H( K, J ) 

Do 64 1=1,12 

Do 64 J=l,12 
Do 64 K=1,12 

64 bE( I ,J)=HYU( ] ,J) /YriY 

Do 65 1 = 1,12 

DO 65 J=l,12 

65 ri ( I , j ) =H { I , J )+AE C I , J )-BE( I , J ) 

DO 66 1=1,12 

GW ( I ) “tva ( J ) 

Gp ( I )=QP1 ( I ) 

66 del i I )^=DEL1 ( I ) 

N = N+1 

IF (N-lOO ) 68 ,6'8 ,57 
68 IF (F2-F1 ) 22,67,22 

67 R=0o25*R 
PR I MT106 » R 

106 FoRMAT( 20X,^THE PRESE’-JT VALUE OF R = «- ,E 1 2 t ) 
KE=KE+1 
PRINT107,KE 
IF(KE-30)23,23,50 
50 RETURN 

End 

$IBFTC OFUNC 

subroutine OFUNC(B,W,G,P) 

DIMENSION B(12) ,G(25) 
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^OMMON/APT/PI ,T ,TEi>FS ,XRl »XR2 9 XR 3 » XRA »F£ND » AN 
v0M,\!0,'\!/r3PT/A, '1 •• A ,;2 9 AN3 » AN A 9 Z? 1 3 ZP2 » ZP3 , 7 pi,. 

PNl=5cO 
PN3i=3oO 
PNA=AoO 

Pl=0o9l2/0o 154 
P2 = 2o0*T/(0ol54*PI-is-^<-2*FEND} 

FB=3o5 

PS=Pl**2/4oO 
A1=B( 1 
A2=B(2)**3 
A3— B( 3 )'^‘^3 
A4=B ( 4 ) **3 

^1= ( B ( 5 ) -lo 0 ) ^^2 

^2 = ( ( 1 ( 6 ) ) / ( B ( 6 ) *"B t 5 ) ) ) ( 5 ) ■“"^^2 

23 = C{B(6)*B(5))*(e(3)<<-B(7)^(loO+8(5))-(B(7)+B(5)))/(B(7)*(B(5)-«'l 
1 )»2.0*B (5)*S(6}+B(6) )-B{5)) 

Z4= ( B ( 8 ) +1 o 0 ) **2/Zl 

w = pS*FB*A 1*Z1*B(9)*-»'2 + PS»FB*A2*Z2*B(9) -if-s-Z+P '3->f A3*FB*Z3-)f-B (11) **2+l 
1FB*A4*24*B( 12 )**2 

W = W/ (10c 0*^e ) 

AN1=(B(5 )-1<.0)*B(9) 

AN1 = ABS( ANI ) 

C22=AN1/(AN1+B(9)-B(6)*B(9) ) 

AN2 = B( 10 )*C2 2*B (6 )*B( 10 ) 

AN2=ABS( AN2 ) 

C 3 3= ( B (7 ) -s-AMl^B ( 9 ) -ANr*( ANl+B (9) ) )/(B(7)*D(<3)*{2 o 0*AN1+B ( 9 ) ) - ( Af 
lB( 9) )**2 ) 

C33=ACs(C33) 

AN3=B(i1)/(C33-1oO) 

AN4=ABS( B (9 )*B{ 12}*(B(8)+1oO) /AMI ) 

2pl=( A n 1-B( 9 ) )/2o0 

2P3=( AN3“6{ 11 ) ) /2o0 
ZP4=(AN4-C(12) )/2oO 
DEN1=B(9 )+ZPl 
DEn3=B(11 )+ZP3 
DEn 4=B ( 12 ) +ZP4 
Sl=(ZPl+2o0)/DEMl 
S3=(ZP3+2 oO) /DEN3 
S4=(ZP4+2oO)/DEN4 
S3=ABS(S3 ) 

S4=ABS(S4) 

Si=ABS£S1 ) 

Si=AMAX1 (S1,S39S4) 

C1=AR5iN{S1) 

C3=ARSIN ( S3 ) 

C4=ARSIN(S4) 

Cl I=pl/Cl 
Cl3=Pl/C3 
Cl4=pl/C4 
D1=ABS(AN1*B(1) ) 

D2:pABS(AN2*B( 2) ) 

D3=ABS( AN3itB( 3) ) 

D4=ABS{AN4#B(4) ) 

DS4=B( 12 )*B(4) 
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DS3 = B(ii)->B(3) 

DS1=B( q ) *B ( 1 ) 

C=0o0033 

X=“0o0009/1000n0 

FST=FEnD* 10000 oO/7oO 
EL = PI*AN/25c.q/12o0 
C^F=:1666o0/0o00i 
DL=2 oO*T*2o2 
Y1=0o154-0c 912/B (q) 

VELSl = DSl4fEL 
W1=DL/D51 
Wl=Wl/pM 

Bl=PI*FE*b( 1 )/25c4 
ERR1=X*VELS1+C 
EFl = CFifERPl 

WDl = 0c05*VELSl*(3iifCFi+Wi) / (0o05*VELSl+S0RT(Bl*CFl+Wl-) )+Wl 
WBl = FST*Bl*YH^PI-!fB(l )/25 o4 
Y3=0o154-0o 912/B ( i 1 ) 

VELS3=DS3*EL 

W3=Dl/DS3 

W3=W3/PN3 

B3 = PI*FB'ifB( 3 ) /25c4 
ERR3 = X*VELS3-hC 
CF3=CF*ERR3 

WD3 = Oo05^^VELS3^:- (B3^^CF3 + W3 ) / ( 0 , 05^-VEL£3+,S?R T ( S’" *CF3+W3 ) )+V/^ 
WB3=FST^-B3*Y3#P I*E(3 ) /25^4 
Y4=0o 154-0o912/B C 12) 

VELS4=DS4*el 

W4=DL/DS4 

W4=W4/pN4 

B4 = pi*Fb^fB ( 4 ) /25 o4 

Err4=x*VElS4+C 

CF4=CF*ERR4 

'^D4=0o 05^VEL3 ' <!■ ( D4*CF4 + W4 ) / ( Oo 05*VELS4 + SQRT ( B4-»-CF4+W4 ) )+W4 

WB4 = FST-J«-B4'»Y4-)'-P 1*3 ( 4 ) / 25c. 4 

G( 1 )=Wpi-l„23*WDl 

G(2)=CI1-3<,0 

G(3) = WB3-lo25i^WD3 

G{4)=CI3-3oO 

G(5)=WB4-1.2:’*WD4 

G (6 ) =C J4-3-I' 

G ( 7 ) =0o85*Dl-D3 
bIFl=ABS ( ( D1-D4 ) /D1 ) 

0IF2=ABS( (01-D2 )/Dl) 

G(8)=0o15-DIFi 
G( 9 ) = 0ul5-DiF2 
G ( 10) =B( 1 )-l oO 
G{11)=R(2)-1 oO 
G(12>=B(3)-1oO 
G(12)=B(4)-1oO 
G( 14) =B( 5 )-4cl95 
6( 15) =4 c.35-B(5) 

G(16)=B(6)-2o3814 
G ( 17) =2c4786-B< 6) 

T G ( 18 ) =B{ 7 )-1o5582 
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G( 19)=1o6218“B(7) 

G(20)=B(8)-5o 8306 
G ( 21 ) =6o 1094-B{ 8 ) 

G(22)=B(9)-15o 99 
G{23)=B{ 10)-15o99 
G(24)=D( 11 )-15o99 
G(25)=B( 12)-16 o0 
P = OuO 

Do 5I=1»25 

5 P=P+lcO/G(I) 

Do 4I=i»25 

IF(G( I ) cGT.OoOGO TO 4 
P=10„0**25 
GO TO 6 
4 CONTINUE 

6 RETURN 

End 

SIBFTC GRAD 

Subroutine grad(B»gw>gpjG) 

Dimension B( 12 ) »g (25) 12 } »gp( 12 ) 

Common/ APT /PI ,T ,TE»FS,XR1>XR2»XP.3 »XR4,Fc:ND,AN 
X=0o0005 

call OFUNC(BiW>G,P) 

X2 = W 
X3 = P 

B ( 1 ) = B( 1 ) + X 

CALL OF'JNC(B,V.'»G,P) 

X4 = W 
X5 = P 

GW(1)=(X4-X2)/X 

GP(1)=(X5-X3)/X 

B(2)=B(2)+X 

B( 1 ) = 6( 1 )-X 

CALL OFUNC ( E ,VSG ,p) 

X6=W 
X7 = P 

GW(2)={X6-X2)/X 
GP(2)=(X7-X3)/X 
6 ( 3 ) = D ( 3 ) +X 
B( 2)=D(2 )-X 
CALL OFUNC ( B >W >G ,P ) 

X31 = W 
X32 = P 

GW(3)~(X31“X2)/X 
Gp (3)=(X32-X3)/X 
B(4)=B(4)+X 
B ( 3 ) =B ( 3 )-X 
Call ofunc( BsWjGsP) 

X41=w 
X42 = P 

GW ( 4) = ( X41-X2 ) /X 
GP ( 4 ) = ( X42-X3 ) /X 
b ( 5 ) = ^ ( S ) + X 
b ( 4) =0 ( 4-)-X 

Call ofui'c(b,w,g.i^‘) 
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X5 1 = W 
X52=P 

GW (5) = (X51-X2)/X 

GP (5) = (X52-X3)/X 

B(6)=B(6)+X 

B(5)=B(51-X 

Gall ofunc(b,w,g,p) 

X61 = w 

X62 = p 

GW (6) = (X61-X2 )/X 
Gp (6) = (X62--X3 )/X 
B( 7)=D(7 )+X 
b (61 =b(6 )-X 

Gall ofuncib.WsGspi 

X 71 =w 

X72 = p 

GW(7)=(X71-X2)/X 

GP (7)=(X72-X3 )/X 

b(8)=D(3 )+X 

B{7)=B(7)-X 

Gall ofunc{B ,w>g,p) 

X81=w 

X82=p 

GW(8)=(XS1-X21/X 

GP(3)=(X82-X3)/X 

B(9)=B(9)+X 

B(8)='3C8)-X 

GALL.0FUNG(B,W,G,P) 

X91=W 

X92 = P 

GW(9)=(X91-X2)/X 
GP{9) = (X92-X3 )/X 
B (9)=B(9 )+X 
B (81=8(0 )-X 

Gall ofunc(B,w,g,pi 

X101=W 

Xl02=p 

GW( 10 1 = (X101-X2 ) /X 
GP( 101=(X101-X2 )/X 

B ( 11) =n( 11 1+x 

B ( 10 1 =B( 10 )-X 

Gall ofunc[B,w»g,p) 

Xlll=w 

Xll2=p 

GW(11)=(X111-X2)/X 
Gp ( 11 ) rr {X112-X3 1 /X 
B ( 12) =B( 12 1+x 
B(111=D( 111-X 

Gall ofunc(B,w»g,p) 

X121=W 

X122=P 

GW( 121=(X121-X2)/X 
GP( 12) = (X122--X3)/X 
RETURN 
EMH) 
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SIBFTC step 

subroutine STEP(B,F,G,F1,DELX»AL 
Dimension bciz) »delx(12)5G(25) ,alp(12) »ald(12) 

COMMON/APT/PI ,T»TE»FS,XRl,XR2»XR3»XR4!.FENDia 

Do 10 1=1,12' 

Y(I)=B{I) 

10 alp( n=-B( n /DELX{ n 

KK=1 

AmiN=ALP(1) 

11 KK=KK+1 
IF(KK»GTo12)G0 to 12 

• IF(AMINoLToALP(KK) )60 to 11 
AMIN=AlP (KK) 

GO To 11 

12 BETA=ABS (AMIN) 
alpha=o<.5*beta ■ 

Cl=ALPHA 

PRINT25»CL 

25 F0RMAT(5X,E10o2) 

IF (CLoGT olaO )CL = 0o5 
22 PF=F 
24 Do 28 1=1,12 
28 G( I )=B{ n+CL*DELX( I ) ■ 

Gall ofunc(C,w,g,p) 

F=w+R*P 

PRINTliO, (C( I ) , 1 = 1,12) ,W,F 
110 Format ( IX, 7( Ei0o4, IX) ) 

IF(PF-F)31,31,32 

31 ■ Cl=0-5*CL 

IF (CL-Oo 0001 *ALPHA) 33,33, 24 

32 Do 34 1=1,12 
3A B<I)=C(I) 

Go TO 22 

33 Call ofunc(B,w,g,p) 

F=W+R«-p' 

Do 40 1 = 1,12 

40 AlD{I)=(B(I)-Y{I) )/DELX( I) 

N = 1 - ■ 

, AMIN=AlD{1) 

41 N=N+1 
IF(NoGTol2)GO TO 42 
IF(AMINoLToALD(N) )GO TO 41 
AmIN=AlD(N) 

GO TO 41 

A2 AlT=ABS( AMIN) 

IF(F1-F)35,35,36 
36 Fi=F 
35 return 

End 

sentry 

2o5 2o5 3oO 3oO 4.28 2o42 1.58 5o96 28.0 28.0 24.0 36.0 




